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FINAL  REPORT 

Summary  and  Analysis 

The  research  reported  here  may  be  divided  into  three  sections.  The  first 
deals  with  the  effect  of  atmospheric  gravity  wave  on  airglow.  The  work  results 
in  a paper  entitled  "On  The  Behaviour  of  Airglow  Under  The  Influence  of  Gravity 
Waves".  (Porter,  Silverman  and  Tuan  (1974)).  The  second  section  deals  with 
"potential"  model  for  gravity  waves  (Yu,  Tai  and  Tuan  (1975)).  Most  of  the  work 
for  this  is  finished.  Upon  completion  we  intend  to  write  a long  paper  on  this  new 
approach  to  gravity  wave  problems.  The  third  section  deals  with  the  non-linear 
effects  of  gravity  waves,  in  particular  their  effects  on  natural  atmospheric 
oscillations  and  airglow,  (Tuan,  Hedinger,  Tai  and  Silverman  (1975)).  A preprint 
is  enclosed  with  this  Final  Report.  We  intend  to  continue  working  on  this  problem 
which  shows  promise  for  handling  a much  more  general  variety  of  non-linear  gravity 
wave  problems. 

I.  On  the  Behaviour  of  Airglow  under  the  Influence  of  Atmospheric  Gravity  Waves 

The  development  of  the  technique  for  treating  the  effect  of  gravity  waves  on 
ionosphere  has  been  given  in  some  detail  in  the  final  report  in  (1973) , (see  also 
Porter  and  Tuan  (1974)).  We  have  since  applied  the  Green's  function  technique 
developed  in  that  paper  to  the  6300  8 01  and  the  5200  8 NI  emissions  in  mid-latitudes. 

Our  procedure  is  based  on  the  assumption  that  these  two  highly  forbidden  emis- 
sion lines  can  only  occur  in  the  F-region  and  the  upper  F-region  where  our  theoretical 
calculations  for  the  behaviour  of  electron  density  profiles  under  the  influence  of 
a fully  ducted  long  period  gravity  wave  (Porter  and  Tuan  (1974))  is  most  valid. 

For  long  period  gravity  waves  the  continuity  equation  for  oxygen  concentration  in 
the  states  reduces  to  the  quasi-equilibrium  equations  of  Petersen,  Van  Zandt 
and  Norton  (1966) . Together  with  the  time-dependent  electron  density  profiles 
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the  emission  rate  for  6300  8 01  may  be  calculated.  A rather  stringent  test  of  our 

theoretical  calculations  has  been  provided  by  the  experimental  observation  of 

Dachs  (1968)  who  has  obtained  not  only  the  6300  8 observations  as  a function  of 

time  but  also  the  simultaneous  f F2  and  h'F2  variations  with  time.  With  a suitable 

o 

choice  of  only  two  parameters,  (the  scale  height  H = 65  Km,  the  scale  height 
factor  p = 1.5),  the  theoretical  results  for  simultaneous  variations  in  6300  8 01, 
fQF2  and  h'F2  appear  to  agree  very  well  with  the  corresponding  experimental  obser- 
vations. Considering  how  few  free  parameters  are  used  in  comparison  with  the 
diversity  of  the  experimental  observations  (both  airglow  and  ionospheric  parameters) , 
the  agreement  must  give  us  confidence  in  our  basic  theory. 

For  the  5200  8 NI  we  can  no  longer  use  the  quasi-equilibrium  theories  of 

2 

Peterson,  Van  Zandt  and  Norton  (1966) , since  the  half  life  of  the  D state  is 

5/2 

very  much  larger  than  the  period  of  the  gravity-wave  oscillation.  This  is  especially 
so  at  high  altitudes  where  quenching  can  no  longer  compensate  for  the  very  slow 
radiative  loss.  By  making  use  of  the  theories  of  Hernandez  and  Turtle  (1969)  which 

neglects  the  diffusion  terms  in  the  continuity  equation  for  atomic  nitrogen  concen- 

2 2 
tration  in  the  D state  but  includes  the  time  dependence,  we  can  calculate  the  D 

nitrogen  state  concentration.  The  effect  of  electrons  on  the  nitrogen  concentra- 
tion is  to  quench  the  optically  forbidden  state  through  superelastic  collisions. 

Thus,  a sinusoidally  varying  electron  density  concentration  produced  by  the  gravity 
waves  would  also  produce  a sinusoidally  varying  5200  8 nitrogen  emission  line. 

Several  interesting  results  have  come  out  of  this  calculation.  (1)  The  impor- 
tance of  electron  quenching  for  this  highly  metastable  nitrogen  state  has  effec- 
tively kept  the  peak  of  the  5200  8 luminosity  profile  relatively  low  even  though 
the  very  long  half  life  may  at  first  seem  to  imply  a high  luminosity  peak.  In  fact 
the  calculated  peak  is  comparable  to  the  6300  8 01  luminosity  peak.  This  is  consis- 
tent with  the  observations  of  Weill  and  Christophe-Glaume  (1967)  who  argued  from 
geometric  considerations  that  the  5200  8 peak  should  only  be  about  12  km  above  the 

__ J 
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6300  X peak  if  the  speed  of  the  travelling  ionospheric  disturbance  is  assumed  to 
be  the  same  for  both  emissions.  (2)  Because  of  the  effect  of  quenching  and  the 
lack  of  0*  and  NO+  molecules  at  high  altitudes  to  produce  airglow,  the  vertical 
movements  of  the  luminosity  profiles  is  considerably  more  restricted  than  the 
movements  of  the  peak  electron  density  profile.  (3)  The  rate  of  change  of  total 
columnar  electron  content  is  always  negative  in  the  absence  of  incident  electron 
flux.  This  means  that  in  the  presence  of  a gravity  wave,  the  oscillation  of  the 
columnar  electron  content  decreases  monotonically  with  height. 

The  result  is  significant  in  that  any  non-monotonic  oscillations  in  the  total 
electron  content  cannot  be  caused  by  gravity  waves.  Thus,  a measurement  of  the 
columnar  electron  content  variation  with  time  can  determine  whether  a gravity  wave 
or  some  incident  electron  flux  is  causing  the  oscillations.  (4)  There  is  always 
a phase  advance  in  the  airglow  fluctuations  with  respect  to  the  change  in  the  total 
electron  content.  The  reason  for  this  is  discussed  in  detail  in  the  enclosed  pre- 
print. 


II.  On  "Potential  Model"  for  Gravity  Waves 
Introduction 

Most  theoretical  studies  on  the  propagation  of  gravity  waves  through  the 
atmosphere  are  classified  into  the  following  categories;  the  ray  theory,  the 
WKB  method  and  the  multi-layer  model.  Of  these  methods  the  multi-layer  full 
wave  approach  (Francis,  S.H.,  (1973),  Klostermeyer , J.,  (1972a, b)  Volland,  M. 


(1969),  Midgeley,  J.E.  and  Liemohn,  H.B.  (1966),  Friedman,  J.P.  (1966))  is  perhaps 
the  most  complete  in  that  the  formalism  allows  the  waves  to  be  scattered  into  dif- 
ferent channels.  For  instance,  a gravity  wave  may  be  scattered  or  transmitted  as 
a heat  conduction  wave  and  vice  versa.  The  basic  assumption  used  for  these  full 
wave  model  calculations  is  that  the  atmosphere  is  horizontally  stratified  and  may 
be  divided  into  thin  isothermal  layers  in  which  we  may  assume  the  various  atmos- 
pheric parameters  to  be  constant.  The  waves  are  then  joined  together  at  the 
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interface  by  appropriate  boundary  conditions.  There  are  a number  of  difficulties 
associated  with  such  an  otherwise  excellent  approach. 

(1)  The  sharp  step  function  change  in  the  parameters  such  as  temperature  etc. 
across  different  layers  would  produce  6-functions  in  their  height  derivatives 
(Hines  (1973) ) . Although  justification  for  some  special  use  of  the  multi-layer 
model  have  been  given  by  a number  of  authors  (e.g.  Pierce,  (1966)),  it  may  be  de- 
sirable to  verify  these  justifications  by  a totally  different  approach. 

(2)  There  has  been  some  uncertainty  about  the  boundary  conditions  to  be  used 
at  the  interface  between  adjacent  layers,  owing  to  the  perturbation  of  the  wave 
on  the  ambient  atmosphere,  (Francis,  1973,  Volland,  1969,  Friedman,  1966). 

(3)  A good  deal  of  uncertainty  seems  to  exist  for  the  so-called  upper  bound- 
dary  conditions.  The  spectrum  for  the  guided  modes  is  dependent  on  whether  this 
upper  boundary  condition  is  assumed  to  be  a free  surface  boundary,  a rigid  sur- 
face boundary  or  the  radiation  condition,  (Friedman  (1966),  Tolstoy  and  Pan  (1970)). 

Basically,  there  is  very  little  one  can  do  about  the  upper  boundary  condition 
as  long  as  one  sticks  to  the  usual  formalism  of  multi-layer  method  as  applied  to 
the  system  of  coupled  first  order  differential  equations.  The  reason  is  that  the 
solution  is  invariably  a plane  progressive  wave  and  in  a numerical  solution  one 
must  follow  such  a wave  up  the  atmosphere  until  an  "upper"  boundary  is  reached. 

In  reducing  the  problem  to  a second  order  Sturm-Liouville  equation  we  may  obtain 
"standing  wave"  solutions  for  the  discrete  eigenstates.  Outside  the  range  of 
the  "potential"  we  may  just  simply  use  the  "outgoing"  boundary  conditions  common^ 
ly  used  in  atomic  and  nuclear  physics. 

(4)  Different  nomenclature  has  been  used  for  waves  considered  not  to  have 
the  same  guiding  mechanism.  A surface  wave,  for  instance,  is  a wave  which  is 
vertically  evanescent  and  is  considered  to  be  guided  by  say  a temperature  dis- 
continuity, (e.g.  the  sharp  temperature  change  at  the  base  of  the  thermosphere) . 
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Such  a ducting  mechanism  has  been  originally  suggested  by  Thome  (1968)  to  account 
for  large  scale  TXD's  in  the  F-region.  The  mechanism  has  been  developed  later 
by  Francis  for  a more  realistic  temperature  profile  (Francis  1973) . Such  a sur- 
face wave  is  not  considered  to  have  the  same  guiding  mechanism  as  for  instance, 
the  S2  mode  at  lower  altitude.  This  mode  is  supposed  to  be  guided  by  two  boun- 
daries instead  of  one,  the  mesopause  and  the  ground  which  form  the  upper  and 
lower  boundaries  of  a wave  guide.  Our  present  approach  should  settle  the  problem 
whether  these  apparently  different  waves  with  different  names  are  in  fact  guided 
by  fundamentally  different  mechanisms. 

(5)  Some  considerable  controversy  exists  for  the  ducting  mechanisms  for 
medium  scale  TID's.  Friedman  (1966)  maintained  that  medium-scale  TID's  are 
caused  by  strongly  ducted  gravity  waves  which  have  rather  low  horizontal  phase 
velocity  (less  than  250  m/sec) . Francis  (1973)  has  maintained  that  gravity  waves 
with  so  low  a horizontal  phase  velocity  cannot  be  ducted  in  the  usual  way.  It 

is  in  fact  refracted  around  the  earth  by  the  changing  direction  of  the  earth's 
radius  of  curvature. 

(6)  So  far  as  we  know  there  are  no  direct  measurements  of  the  presence  of 
gravity  waves  at  F-region  heights.  Usually,  the  waves  are  detected  indirectly 
through  its  effect  on  the  electron  density  concentration  which  can  then  be  mea- 
sured by  means  such  as  the  incoherent  backscatter.  A knowledge  of  the  phase  front 
of  such  waves  can  help  us  determine  the  degree  of  ducting.  In  general,  it  is  not 
easy  to  deduce  the  direction  of  the  phase  front  from  the  phase  of  the  TID.  Porter 
and  Tuan  (1974)  have  shown  that  if  a purely  vertical  phase  front  is  assumed  for 
the  gravity  wave,  the  TID  would  have  a nearly  horizontal  phase  front  in  the  F- 

I 

region  peak.  In  fact  phase  variation  with  height  is  almost  identical  to  Fig.  14 
in  Francis'  paper  (Francis  1973)  (see  also  Figure  7,  Porter  and  Tuan  1974).  Thus, 
just  because  the  phase  front  of  TID  may  be  nearly  horizontal  at  the  F-region  peak. 


does  not  imply  that  the  gravity  wave  has  a horinzontal  phase  front.  Thome  (1968) 
has  used  a step  function  temperature  profile  to  obtain  a single  fully  ducted  mode. 
Francis  (1973)  has  used  the  1965  Cospar  Model  and  obtained  three  partially  ducted 
modes  (for  long  period  gravity  waves) . 

In  view  of  all  these  difficulties  some  of  which  apparently  arise  from  the 
multi-layer  model  itself,  we  will  investigate  the  behaviour  of  atmospheric  gravity 
waves  by  a totally  different  approach.  We  will  treat  the  gravity  wave  propagation 
through  the  atmosphere  as  a multi-channel  scattering  problem  with  effective 
"potentials"  for  each  channel.  For  this  report  we  will  restrict  ourselves  to 
altitudes  below  500  Km  where  we  may  neglect  viscosity.  The  gravity  waves  may  then 
be  scattered  into  the  heat  conduction  wave  channel  and  vice  versa.  The  problem  be- 
comes very  similar  to  scattering  in  atomic  or  nuclear  physics.  In  first  approxi- 
mation we  will  decouple  the  heat  conduction  wave  channel  from  the  gravity  wave 
channel  by  assuming  a lossless  atmosphere.  For  the  gravity  wave  channel,  we  obtained 
the  "potentials"  for  Summer  and  Winter  Seasons  as  well  as  different  exospheric 
temperature  from  the  1972  Cospar  Models.  The  problem  may  then  be  treated  as  a 
scattering  process  (free  modes) , a resonance  process  (partially  guided  modes)  or  a 
bound-state  problem  (fully  ducted  modes) . 

We  will  investigate  the  behaviour  of  this  potential  for  different  gravity  wave 
frequencies  as  well  as  for  different  seasons  and  different  exospheric  temperatures. 

We  will  also  determine  the  number  of  bound  (fully  guided)  and  resonance  (partially 
guided)  states  together  with  the  behaviour  of  the  pressure  variation  as  a function 
of  height  for  those  states.  We  will  also  show  that  for  long  period  waves  the 
boundary  conditions  together  with  the  potential  imply  the  existence  of  a lowest 
ground  state  belonging  to  the  lowest  eigenvalue.  Since  the  eigenvalue  is  in  units 
of  inverse  horizontal  phase  velocity,  there  exists  an  upper  limit  to  the  horizontal 
phase  velocity  for  guided  modes.  The  free  modes  also  have  an  upper  limit  (somewhat 
higher)  which  depends  on  the  exospheric  temperature.  We  will  also  develop  an 
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approximate  semi-empirical  dispersion  relation  for  a lossless  atmosphere  for  alti- 
tudes higher  than  130  Km  where  the  "potential"  becomes  very  flat  and  structureless. 
With  this  dispersion  relation  we  hope  to  show  that  an  F-region  gravity  wave  can  to 
some  extent  be  ducted  by  the  curvature  of  the  earth  and  that,  in  fact,  the  sharp 
temperature  change  at  the  base  of  the  thermosphere  as  given  by  the  Cospar  models 
cannot  produce  a sufficiently  deep  "Potential  well"  at  the  F-region  to  guide  a 
gravity  wave. 

We  will  consider  just  how  sharp  the  temperature  change  will  have  to  be  in 
order  to  produce  a well  deep  enough  to  guide  the  F-region  gravity  wave. 

Theory 

As  a first  step,  we  shall  consider  only  two  channels,  the  gravity  wave  and 
the  heat  conduction  channel.  This  is  not  a bad  approximation  (Volland,  (1969))  for 
altitudes  below  500  Km.  It  can  be  shown  that  for  horizontally  stratified  atmosphere 
the  coupled  hydrodynamic  equations  can  be  reduced  to  the  following  coupled  second 
order  eigenvalue  equations : 

.LJM.  + a liM.  + 3 (Ap)  - K2  (Ap)  = A (AT)  + B --(-A-Tj-  (1) 

3z  3z  X 3z 


_3_ 

3z 


* 3 (AT) 
3z 


+ D ( AT) 


k K2  (AT)  = y (Ap)  + 6 

X d Z 


where  a,8,Y/<5  ,A,B,D  are  known  functions  of  temperature,  mean  molecular  mass,  ratio 
of  specific  heats  and  thermal  conductivity.  K^  is  the  horizontal  wave  number  and 
K is  the  thermal  conductivity.  Equation  (1)  contains  both  the  gravity  wave  solution 
and  the  heat  conduction  wave  solution.  Each  type  of  wave  would  have  two  linearly 
independent  solutions  corresponding  to  upward  and  downward  propagation.  In  the 
case  when  k = 0 we  no  longer  have  heat  conduction.  AT  may  then  be  immediately 
eliminated  and  we  would  obtain  a single  second  order  eigenvalue  equation  whose 
solutions  are  just  the  gravity  waves. 
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Since  the  boundary  conditions  for  the  dependent  field  variables  (Ap,AT,eec. ) 

are  just  the  continuity  of  the  field  variable  together  with  their  hydrodynamic 

derivatives,  (Thome,  (1968)),  they  are  essentially  the  same  as  the  quantum  mechani- 

2 

cal  boundary  condition.  Thus,  the  eigenvalues  which  are  just  the  horizontal 
wave  numbers  may  form  a discrete  set  corresponding  to  "ducted"  modes  or  a continuum 
corresponding  to  "free"  modes. 

( .1 ) Discrete  (bound  or  resonance)  States 

When  the  thermal  conductivity  k is  set  equal  to  zero,  equation  (1)  will  immedi- 
ately reduce  to  the  following  equation 


c 


dz  2(z)  dz 


v ( z ) + X' 


U = 0 


(2) 


where  y = to 


exp[-  y jf'/f  dzJAp 


f y , 2 2. 

f = - (<ob  - to  ) 


u.'b  = Brunt-Vaisala  frequency 
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2 2 
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\ C / “ 0) 
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Immediately  from  equation(2),  we  learn  that  (1)  for  ducted  or  paitially  ducted 

2 1 

modes,  since  \p  obeys  the  usual  boundary  conditions,  the  eigenvalue  X^  = ^2  will 

phx* 

be  discrete  which  means  that  the  horizontal  phase  velocity  will  also  be  discrete. 
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Furthermore,  there  is  clearly  a "ground"  state  corresponding  to  maximum  horizontal 
phase  velocity. 

2 2 

(2)  For  long  period  gravity  waves,  >>  to  , the  "potential"  becomes  indepen- 

dent of  the  frequency  of  the  gravity  wave.  In  practice  if  the  period  of  the  gravity 
wave  is  greater  than  30  minutes,  the  "potential"  is,  for  all  practical  purposes, 
independent  of  u).  This  means  that  for  such  waves  there  is  no  longer  significant 
dispersion.  It  is  easy  to  show  that  for  an  isothermal  atmosphere,  and  constant 
sound  speed  profile  equation  (2)  would  reduce  to  the  Hines'  dispersion  relation  if 
we  substitute  in  a plane  wave  solution.  The  profiles  of  temperature,  mean  molecular 

weight,  ratio  of  specific  heats  and  gravitational  constant  are  taken  from  the  1972 

* 

Cospar  models,  and  the  U.S.  stand  and  atmosphere  (1965).  Fig.  1 shows  the  tempera- 
ture profiles  for  the  Summer  Season  with  exospheric  temperature  ranging  from  600°K 
to  2200°K.  Fig.  2,  3 and  4 show  the  "potentials"  for  long  period  gravity  wa^es 
(>  30  min.)  when  they  are  no  longer  dependent  on  gravity  wave  frequency.  They 
correspond  to  exospheric  temperatures  of  600°K,  1800°K  and  2200°K  respectively. 

The  dashed  curve  is  the  Summer  "potential"  while  the  solid  line  is  the  Winter 
"potential" . 

These  potential  curves  above  can  already  provide  us  with  a lot  of  information. 
One  at  once  sees  that  there  are  no  strictly  fully  ducted  modes  since  clearly  only 
"resonance"  states  can  exist  for  such  potentials.  For  long  period  gravity  waves, 
there  is  clearly  an  upper  limit  to  the  horizontal  phase  velocity  as  well  as  an 
upper  limit  to  the  horizontal  phase  velocity  of  the  guided  modes.  We  should 
emphasize  here  that  in  neglecting  viscosity  and  heat  conduction,  our  "potential" 
models  are  good  approximations  only  below  150  Km.  Above  150  Km,  the  dissipative 
losses  become  important.  It  is  possible  to  show  that  such  losses  form  an  imaginary 
part  in  addition  to  our  real  "potential"  given  by  equation  (3) . We  may  then  in 
general  use  the  optical  model  treatment  for  the  entire  problem. 


12 


f 


We  will  for  the  moment  consider  only  the  real  "potential"  given  by  equation 
(3) . If  we  use  a Winter  temperature  profile  with  an  exospheric  temperature  of 
1500°K,  the  potential  for  a gravity  wave  with  a two-hour  period  is  given  by  Fig.  5. 
The  discontinuity  at  the  ground  surface  provides  a 6-function  "potential"  which 
may  be  replaced  by  the  boundary  condition. 


ii  | = (l  ( 12.  \ . £'  A 

8z  lz=0  LAcA  c c2j 


(4) 


z=o 


The  important  feature  of  our  "potential"  is  that  it  becomes  quite  flat  and 
constant  above  200  Km.  The  solution  above  this  height  can  be  approximated  very 
well  by  plane  waves.  The  problem  of  determining  ducted  or  partially  ducted  modes 
thus  reduces  to  determining  the  bound  or  resonance  states  for  atomic  and  nuclear 
physics.  We  may  hence  use  phase  shift  analysis  to  determine  the  resonance  states. 
The  phase  shift  may  be  computed  from  the  usual  equation 


K cot  (K  R + 6 (X2)  ) = ^ 
z z \jj 


(5) 


R 


where  R is  any  altitude  in  the  region  where  the  potential  is  effectively  flat.  A 


resonance  occurs  whenever  6 goes  through  — . Fig.  6 shows  a plot  of  the  phase 


shift  computed  from  equation  (5) . The  curve  shows  that  the  phase  goes  through  a 


2 2 

very  sharp  change  of  tt  radians  at  X = 10.43  and  X =14.33  corresponding  to  horizon- 


tal phase  velocities  of  309.7  and  264.2  meters  per  second  respectively.  This  agrees 


very  closely  with  the  and  modes  of  Francis  (1973) . A rather  faint  suggestion 


of  a resonance  occurs  for  X ^26  but  the  resonance  is  too  broad  to  be  considered 


as  anything  much  more  than  background.  Something  similar  has  been  discussed  by 

2 


Francis  with  the  multi-layer  method.  The  very  sharp  resonance  at  X =10.43  is  the 
well  known  Lamb  mode.  As  can  be  seen  from  the  location  of  the  two  resonances  in 
Fig.  5 the  Lamb  mode  is  not  really  fully  ducted.  In  fact,  there  are  no  fully  ducted 
modes  as  already  mentioned  not  even  the  Lamb  mode  which  has  the  lowest  eigenvalue. 
Considering  the  fact  that  we  have  not  used  similar  temperature  profiles  as  Francis 
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and  that  we  are  only  considering  the  real  part  of  the  potential  the  agreement  would 
be  another  confirmation  that  the  imaginary  part  of  the  potential  (the  dissipation) 
is  relatively  unimportant  below  150  Km.  We  will  hence  identify  our  resonance  as 
the  and  S 2 modes.  The  eigenfunction  for  the  mode  is  given  in  Fig.  7 which 
shows  the  very  high  amplitude  near  the  ground  level.  Above  150  Km  when  the  poten- 
tial is  very  flat  with  no  structure  whatsoever,  the  amplitude  of  the  eigenfunction 
becomes  exceedingly  small.  However,  there  is  clearly  a small  finite  possibility 
for  leakage  and  the  eigenfunction  is  definitely  that  of  a resonance  state.  Fig.  8 
shows  a plot  of  the  S 2 mode  eigenfunction  \p2-  Once  again,  the  amplitude  becomes 
smaller  above  150  Km  and  is  quite  large  below  that  height  level.  However,  there  is 
clearly  considerably  more  leakage.  \p2  has  large  amplitudes  at  about  50  Km  where 
there  is  a deep  well  (Fig.  5)  formed  by  the  hump  at  the  stratopause.  However  the 
most  important  guiding  mechanism  for  this  wave  still  comes  from  the  6 -function  po- 
tential at  the  ground  level.  For  this  reason,  the  partially  guided  gravity  waves 
for  both  the  and  S 2 modes  should  be  readily  observable  at  the  ground  level. 

(4)  The  Guided  Mode  at  the  F-Region. 

At  F-region  altitudes  dissipative  losses  due  to  heat  conduction,  viscosity  and 
ions  drag  become  very  important.  Perhaps  the  most  important  of  these  losses  is 
heat  conduction  ,Volland  (1969),  Francis  (1973a)). 

Since  at  the  moment  we  are  only  concerned  with  the  real  part  of  the  potential, 
we  will  consider  the  effect  of  the  temperature  structure  as  well  as  the  variation 
in  the  mean  molecular  mass  on  long  period  gravity  waves  with  high  horizontal  phase 
velocity. 

There  is  a general  belief  that  the  sharp  temperature  drop  at  the  base  of  the 
thermosphere  can  be  considered  as  some  form  of  discontinuity  which  can  guide  a 
surface  wave.  Thome  (1969)  has  used  a step  function  temperature  profile  to  obtain 
a single  fully  guided  mode.  Later  Francis  (1973)  by  using  the  1965  Cospar  model 
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together  with  assumptions  on  losses  due  to  heat  conduction  and  viscosity  succeeded 
in  obtaining  three  partially  ducted  modes.  Both  Thome  and  Francis  used  the  multi- 
layer method,  although  in  the  case  of  Thome  there  were  only  two  layers,  since  a 
step  function  temperature  profile  was  used. 

Within  the  physical  framework  of  a potential  model,  a step  function  temperature 
profile  is  just  a 6-function  potential.  It  is  well  known  that  a 6-function  poten- 
tial can  have  a single  bound  state  with  finite  eigenvalue.  Thus,  the  single  mode 
of  Thome  is  not  surprising.  Furthermore,  a well  behaved  solution  for  a 6-function 
potential  must  consist  of  decaying  exponentials  which  implies  that  the  wave  numbers 
are  purely  imaginary. 

Since  this  is  just  the  definition  of  a surface  wave  ducted  by  a single  discon- 
tinuity, we  see  that  within  the  present  physical  framework,  there  is  no  need  to 
consider  a surface  wave  as  something  separate  from  other  guided  waves. 

A glance  at  Fig.  2,  3 and  4 immediately  reveals  the  fact  that  the  F-region  well 
if  it  exists  is  very  shallow  indeed.  Part  of  the  reason  why  it  is  so  shallow  is 
that  the  "potential"  as  given  by  equation  (3)  really  depends  on  derivatives  of  the 
speed  of  sound  c rather  then  the  temperature  profile.  The  variation  of  the  mean 
molecular  mass  with  altitude  is  such  that  the  sound  speed  profile  has  a much  more 
gradual  slope  than  the  temperature  profile. 

However,  for  high  exospheric  temperatures,  it  is  possible  to  show  that  an  ex- 
ceptionally shallow  but  broad  well  does  exist  at  the  F-region  which  corresponds  to 
the  temperature  rise  above  the  mesopause.  Such  a potential  well  can  be  very  accurately 
approximated  by  the  Morse  Potential  given  by: 


V = VQ  + a [e_2a (z-zo)  - 2e-°(z-zo>] 

inhere  zq  = height  at  which  the  minimum  of  the  well  occurs 

V = some  scaling  constant. 
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The  quantities  a and  a are  the  two  Morse  potential  parameters.  One  can  then  show 

that  the  criteria  for  the  existence  of  bound  states  in  terms  of  the  Brunt  frequency 

is  given  by: 

J , 2 2. 

va  (wb  - a)  ) ^ 

N = > •=• 

a 2 

where  u)  is  the  frequency  of  the  gravity  wave.  Table  I shows  the  values  for  N for 
long  period  (>  -j  hour)  gravity  waves  and  exospheric  temperatures  of  2000°K  and 
2200°K  respectively. 

TABLE  I 


Exospheric  Temp. 

2000° 

2200° 

Z 

o 

272  Km 

278  Km 

I 

a 

1/70  Km 

1/74  Km"1 

A 

7.396  x 10_2Km_2sec"2 

8.39  x 10  2Km  2sec  2 

V 

o 

6.86x10  1Km  2sec2 

5.957  x 10  ^'Km  2sec2 

I 

N 

0.133 

0.146 

It  is  clear  that  both  values  of  N are  less  than  0.5  the  minimum  required 

for  a bound  state.  For  lower  exospheric  temperatures  N becomes  smaller  still. 

In  order  to  determine  just  how  sharp  the  temperature  drop  will  have  to  be  at 

the  base  of  the  thermosphere  in  order  to  produce  a well  sufficiently  deep  to  guide 

a gravity  wave  we  used  an  analytic  temperature  profile  given  by  Fig.  9a.  Fig.  9b 

shows  the  corresponding  well  with  a minimum  at  about  140  Km.  It  is  clear  that 

such  a well  can  only  produce  fully  ducted  modes  with  horizontal  phase  velocities 

2 

greater  than  1000  meters/sec  (X  < 1).  We  are  not  aware  of  observations  which 
support  this. 

To  explain  Thome's  (1969)  observations  which  gave  pretty  convincing  evidence 
that  at  least  a partially  ducted  wave  at  the  F-region  altitudes  exist  with  horizontal 
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phase  velocities  of  the  order  of  700  m/sec  and  a period  of  about  2 hours,  we  must 


point  out  that  the  real  part  of  the  "potential"  only  becomes  predominant  below 


150  Km.  Above  200  Km  the  imaginary  part  of  the  potential  becomes  much  more  impor- 


tant than  the  real  part.  If  we  confine  our  attention  for  the  time  being  to  below 


150  Km,  we  see  immediately  from  the  potential  that,  for  X <8  which  corresponds 


to  v ^ > 350  m/sec,  there  is  little  gravity  wave  below  120  Km.  Fig.  10  shows  a 


diagram  of  a gravity  wave  with  a 2 hour  period  and  a horizontal  phase  velocity  of 


707  m/sec.  which  is  rather  similar  to  Thome’s  (1968)  observed  gravity  waves.  We 


see  that  there  is  relatively  little  wave  amplitude  below  130  Km.  This  is  in  good 


agreement  with  Thome's  experimental  observation  of  the  gravity  wave  on  May  13-14 


and  on  June  9-10,  1964.  The  May  13-14  wave  has  horizontal  phase  velocity  of  700  m/s« 


and  induces  imperceptible  oscillation  at  134  Km  where  the  electron  density  is  still 


relatively  significant.  The  June  9-10  wave  has  a horizontal  phase  velocity  of 


600  m/sec  during  the  day  and  again  induces  no  oscillation  at  131  Km.  From  Fig.  10, 


we  see  that  waves  with  such  horizontal  phase  velocities  do  not  have  very  much  ampli- 


tude below  130  Km.  In  fact,  the  gravity  wave  becomes  evanescent  below  120  Km. 


These  results  are  also  in  good  qualitative  agreement  with  the  Gq ' , ' , and  G^' 


pseudomodes  of  Francis  (1973) 


Above  150  Km  we  can  no  longer  ignore  the  imaginary  part  of  the  potential.  To 


a first  approximation,  we  may  use  the  WKB  approximations  which  can  immediately  show 


that  the  imaginary  part  produces  severe  amplitude  damping  with  height  as  well  as  a 


change  of  phase  in  the  eigenfunction  4i.  The  eigenvalues  X of  the  eigenstates 


will  also  change.  However,  since  imaginary  part  W of  our  potential  is  so  large 


at  higher  altitudes  that  we  can  no  longer  ignore  second  order  terms  in  the  scatter- 


ing amplitude  A.  We  may  write 


A = ~ vtGovt  + ••••IV 


where 


A 
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= v + iw, 

and  Gq  is  the  Green's  function  obeying  the  same  boundary  condition  as  ip  at  the 
origin  and  an  outgoing  boundary  condition  at  z + » , Gq  satisfies  the  equation 

Q<Z)  dl  + X2  GlZ'Z'’X2)  = -5(2-2')  . 

The  second  order  term  contains  four  terms,  one  of  which  is  the  real  term 
WGqW  which  can  be  very  important  and  may  produce  resonance  states.  Since  W -*■  0 
as  k + 0,  we  know  that  such  resonance  states  corresponding  to  partially  guided 
modes  are  strictly  produced  by  dissipation. 

In  summary  for  the  F-region  long  period  gravity  wave,  we  are  in  agreement 
with  both  the  experimental  observations  of  Thome  as  well  as  the  theoretical  treat- 
ment of  Francis  below  140  Km.  Above  this  height  level  we  are  not  in  any  real  dis- 
agreement with  Thome  since  we  would  have  obtained  a single  fully  ducted  mode  if 
we  had  used  a step  function  temperature  profile.  However,  we  are  in  sharp  dis- 
agreement with  Francis  (1973)  who  has  obtained  three  partially  ducted  modes  which 
he  has  labelled  G ' , G ' , G ' and  has  furthermore  maintained  that  they  are  essen- 
tially  produced  by  the  temperature  discontinuity  at  the  base  of  the  thermosphere. 
Since  we  have  both  used  realistic  temperatures  as  well  as  speed  of  sound  profiles, 
it  is  difficult  to  see  how  such  differences  may  be  easily  reconciled.  As  already 
mentioned,  the  only  way  we  can  obtain  ducted  modes  is  through  the  second  order 
effect  of  the  imaginary  part  of  our  potential.  Since  W,  the  imaginary  part  of 
the  potential,  is  large,  it  is  certainly  possible  that  some  kind  of  ducting  cor- 
responding to  a resonance  state  may  be  feasible.  However  such  resonance  states 
are  not  caused  by  the  temperature  change  at  the  base  of  the  thermosphere. 

(5)  The  effect  of  the  earth's  curvature  on  F-region  gravity  waves. 

Although  the  temperature  gradient  appears  to  be  insufficient  for  ducting 


longer  period  atmospheric  gravity  waves  in  the  F-region,  there  are  means  other 


r 


1 


18 

than  atmospheric  structure  which  may  be  able  to  produce  sufficient  ducting  to  ex- 
plain the  observed  behaviour.  Practically  all  observed  F- region  gravity  waves 
are  ducted  only  about  1/4  of  the  way  around  the  earth.  The  distance  is  of  the 
order  of  10,000  Km.  Fortunately,  for  the  F-region,  the  "potential"  is  very  nearly 
flat  and  cd^  also  becomes  a very  slowly  varying  function  of  height.  Thus,  we  may 
to  a first  approximation  assume  that  Q( z)  and  v(z,w)  are  both  constant.  From 
equation  (2)  we  may  immediately  derive  the  following  semi-empirical  dispersion 


relation  for  the  F-region 

kz2  = kx2(  -J-  _1)  ‘ 
0) 


(,  2 - ID2)  V 
b 


(6) 


where  V is  a constant  and  may  be  taken  from  the  flat  part  of  the  "potential"  curves. 

We  may  now  adopt  a similar  procedure  as  Francis  (1972) . It  is  possible  to  show  that 

the  group  velocities  are  given  by: 
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If  we  introduce  x as  a variable  defined  by 
x = 1^0 

where  Rg  is  the  radius  of  earth  and  0 is  the  angle  with  respect  to  say  the  axis 

of  rotation,  and  z as  the  vertical  height  normal  to  the  earth's  surface,  then 
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dz 

dx 


v 

gz  _ 


gx 


u>2 

kxHr  - 15 

CD 


(8) 


Since  V is  a constant  in  equation  (6) , we  may  immediately  obtain  the  follow- 


ing equations  for  Kx  and  K^ 
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dK  K 

z x 


dx 


*E 


dK  K 
x z 


dx  re2  wb  ” 10 

The  above  equations  are  identical  to  those  of  Francis  (1972) . 


Solving  for  and  Kz»  we  obtain 


[' 


0) 


-]  — -[-  “ 


. 2 2,1/2  J r_  ‘ 2 2, 1/2J  R 

(w,  - a)  ) E (a),  - a)  ) 

D D 

Kx  = A e + D e 


E 


l- 


b) 


/ 2 , 2 V.2_  ...2,1/2 


■]  — ~ 


0) 


w (U)^-  rE  (U)2-  w2)1/2  *E 

K = — ^ (A  e b - D e b ) 

Z b) 


(9) 


(10) 


Since  D > A and  the  exponent  is  much  less  than  one  for  x 10,000  Km,  is  negative. 
There  is  hence  always  a strong  forward  tilt  as  shown  in  Fig.  11.  If  we  assume 
that  the  ray  path  begins  at  160  Km,  it  will  reach  a height  of  only  about  310  Km 
after  travelling  10,000  Km  i.e.  1/4  of  the  way  around  the  earth. 

The  gravity  wave  here  is  assumed  to  have  a two  hour  period  and  a horizontal 
phase  velocity  of  723  m/sec.  Since  these  numbers  are  similar  to  the  waves  observed 
by  Thome  (1968) , it  would  seem  not  beyond  the  realm  of  possibility  that  the  earth's 
curvature  can  be  an  important  contributing  factor  to  F-region  gravity  wave  ducting. 


(6)  "Potential"  Model  and  the  Hines  Dispersion  Relation. 

In  regions  where  the  Brunt-Vaisala  frequency  is  a slowly  varying  function  of 
height,  the  "potential"  model  reduces  to  the  Hine's  theory  for  small  horizontal 
phase  velocity;  (large  X 2).  From  equation  (2)  we  see  that  if  X 2 » |v|,  (2) 

X x 

reduces  to 


d + , b ,,  „ 2 


dz 


( -J~u  V* = 0 


(12) 
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Substitution  of  plane  wave  solution  e z in  (12)  will  immediately  yield  the  Hines 
dispersion  relation. 

-2  2 

Generally  speaking,  the  "potential"  is  of  the  order  of  25  Km  sec  which 
corresponds  to  a horizontal  phase  velocity  of  200  m/sec.  This  means  that  for 
Vphx  <<  200  m/sec,  the  Hines  model  would  be  rather  good  in  regions  where  wB  is  a 
slowly  varying  function  of  height. 

To  compare  the  various  solutions  we  have  shown  a plot  of  the  and  S2  modes 
together  with  the  "free"  gravity  wave  solution  for  horizontal  phase  velocities  of 
707.1  m/sec  and  213.2  m/sec  in  Fig.  12.  As  can  be  seen,  it  is  difficult  to  plot 
them  all  on  the  same  scale.  The  "free"  gravity  wave  corresponding  to  213.2  m/sec 
is  already  approaching  a simple  sine  wave.  It  is  interesting  to  note  that  in  the 
"potential"  picture  one  sees  immediately  that  all  gravity  waves  with  horizontal 
phase  velocity  higher  than  310  m/sec  must  necessarily  be  F-region  waves.  They  must 
all  become  evanescent  below  about  120  Km.  The  707.1  m/s ec  gravity  wave  clearly 
illustrates  this. 

We  have  so  far  confined  ourselves  to  investigating  gravity  waves  with  periods 
larger  than  30  min.  when  the  "potentials"  do  not  vary  significantly  with  frequency. 
For  gravity  waves  with  shorter  periods  than  30  min.  the  "potential"  curves  are 
drastically  altered.  Fig.  13  provides  an  example  for  a gravity  wave  with  only 
a 7.5  min.  period.  The  investigation  of  such  types  of  waves  is  beyond  the  scope 
of  the  present  report. 


Ill . On  Non-Linear  Effects  of  Gravity  Waves 

There  are  many  airglow  oscillations  showing  a small  period  (8-14  minutes) 
oscillation  superimposed  over  a long  period  oscillation  of  the  order  of  hours. 
Most  of  the  data  available  are  given  by  the  5577  8 01  airglow  emission  and  to  a 
lesser  extent  the  6300  8 red  line  emission.  Sample  examples  of  such  data  will  be 
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to  state  the  principal  objectives  of  our  proposal.  The  airglow  data  may  be  summar- 
ized as  follows. 

(1)  Most  of  the  data  are  not  tied  to  any  known  impulsive  energy  sources 
such  as  an  explosion. 

(2)  Some  of  the  data  shows  that  the  small  period  oscillation  follows 
the  long  period  oscillation  in  time  sequence  and  is  in  fact  slowly  excited  as  the 
long  period  oscillation  increases  in  intensity. 

(3)  There  appears  to  be  some  variation  in  the  period  of  the  small  period 
oscillation  with  the  undulation  of  the  long  period  oscillations. 

There  are  theoretical  analyses  which  may  possibly  be  invoked  to  explain  these 
results  (although  the  original  authors  may  not  quite  have  had  such  a situation  in 
mind  when  they  proposed  their  theory) , for  instance  the  treatments  of  Row  (1967) 
and  Liu  and  Yeh  (1971).  Row's  analysis  deals  with  point  impulsive  sources.  Liu 
and  Yeh's  treatment  deals  with  a generalization  of  the  point  source  to  include  spa- 
tial and  temporal  extensions.  They  were  also  able  to  show  that  at  far  distances 
("v  1000  Km)  away  from  the  source  we  can  obtain  a short  period  Accoustic  frequency 
oscillation  superimposed  over  the  long-period  gravity-wave  oscillation.  While  we 
have  no  doubt  that  their  treatment  is  sound  we  do  have  some  reservations  over  whether 
it  can  be  used  for  explaining  the  observed  results  for  the  following  reasons : 

(1)  These  theoretical  treatments  do  require  some  kind  of  impulsive  energy 
sources. 

(2)  The  results  of  Liu  and  Yeh  (1971)  show  that  at  far  distances  from  the 
source  (^  1000  Km)  the  short-period  (accoustic  frequency)  wave  arrives  first  fol- 
lowed by  the  long  period  gravity  wave. 

(3)  In  both  these  treatments  an  inviscid  lossless  atmosphere  is  assumed.  If 
we  were  to  use  such  a theory  to  explain  the  airglow  oscillation  we  must  hence  assume 
that  both  the  gravity  wave  and  the  accoustic  frequency  wave  originate  from  far  dis- 
tances. 
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(4)  The  observed  short  period  oscillation  occurs  only  when  the  long  period 
oscillation  has  a large  amplitude.  Liu  and  Yeh's  (1971)  treatment  is  essentially 
a linear  treatment  and  does  not  allow  a sharp  cut-off  of  the  short  period  oscilla- 
tion at  certain  amplitudes  in  the  long  period  oscillation. 

We  believe  that  the  observed  airglow  oscillations  just  mentioned  are  typical 
examples  of  non-linear  coupling  between  a gravity  wave  and  the  local  atmosphere. 

We  propose  the  following  simple  model  which  we  intend  to  later  develop  into  a full 
scale  theory  based  on  hydrodynamic  equations. 

Generally  speaking,  the  hydrodynamic  equations  are  coupled  non-linear  partial 
differential  equations  for  dependent  field  variables  such  as  pressure  p(x,t) 
density  p (x  ,t)  and  velocity  v(x,t)  . From  a mathematical  point  of  view  it  is 
very  difficult  to  handle  non-linear  partial  differential  equations  and  from  a 
physical  point  of  view  it  is  also  a lot  easier  to  visualize  a lump  of  fluid  in 
oscillation  rather  than  dealing  with  field  variables.  Our  first  step  is  therefore 
to  go  from  a field  picture  into  a particle  (lumped  mass)  picture. 

In  general,  if  f = f(z,t)  is  any  field  observable  (f  may  be  the  pressure, 
density  or  temperature),  f(z(t,zQ),t)  is  the  same  observable  of  a particular  element 
of  fluid  (a  particle)  identified  by  its  initial  position  zQ  at  t = 0.  Here  z = 


z(t,zQ)  is  of  course  the  trajectory  of  the  fluid  element.  It  is  possible  to  show 


that 


3f  (z,t) 
3t 


or 


3f  (z,t) 
3z 


may  also  go  into  their  corresponding  particle  picture  if 


we  substitute  z = z(t,zQ)  in 


3f  (z,t) 
3t 


or 


3f  (z,t) 
3z 


Furthermore,  it  makes  no  differ- 


ence whether  we  substitute  in  z(t,zQ)  first  and  then  differentiate  or  differentiate 
first  and  then  substitute.  In  the  former  case  ^ ■ is  just  the  change 

in  f with  respect  to  a change  in  the  particle  trajectory  while  keeping  the  time 
fixed.  The  same  applies  to  higher  order  derivatives. 

We  shall  for  the  present  report  neglect  heat  conduction  and  only  use  phenomeno- 
logy for  viscosity.  We  shall  also  assume  that  variations  along  the  horizontal 
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directions  may  be  neglected.  It  is  then  easy  to  show  that  the  hydrodynamic 
derivative 


Df 

Df 


'z  = z (t,  zq) 


3f  3f 

3t  + 2 3z 


We  will  use  the  above  procedure  to  convert  Euler's  field  equation  into  a 
particle  equation.  As  a first  step,  we  will  derive  the  equation  for  the  natural 
oscillation  of  the  atmosphere  in  an  unperturbed  hydrostatic  background.  The 
usual  method  used  for  such  a derivation  has  always  been  within  the  framework  of 
a particle  picture  in  which  the  "particle"  is  an  element  of  fluid  contained  in  a 
weightless  invisible  plastic  bag.  In  order  to  later  incorporate  the  gravity  wave 
field  as  a background  in  addition  to  the  hydrostatic  background  and  consider  the 
coupling  between  the  gravity  wave  and  the  Brunt  oscillation,  it  is  much  more  natural 
to  derive  the  Brunt  oscillation  from  a field  picture. 

The  field  equations  can  then  be  converted  into  "particle"  equation  which  is 
just  the  equation  for  simple  harmonic  motion.  We  can  then  use  very  similar  pro- 
cedure to  derive  the  Brunt  oscillation  with  the  gravity  wave  field  as  background. 
After  conversion  into  a "particle”  equation  we  would  have  a Hill's  type  of  equation. 

We  will  show  that  the  time  dependent  parameter  will  have  a fundamental  fre- 
quency equal  to  the  gravity  wave  frequency  and  that  parametric  excitation  of  the 
natural  frequencies  can  occur.  Since  the  velocity  gradient  term  in  the  hydro- 
dynamic  equation  is  generally  considered  the  most  important  non-linear  term  we 
will  incorporate  this  into  our  original  Euler's  field  equation.  The  term  becomes 
a driving  term  for  the  "particle"  equation  (10)  in  the  enclosed  preprint. 


Theory 

We  begin  by  considering  the  displacement  of  a fluid  element  from  its  equili- 
brium position  zq  in  a hydrostatic  background  atmosphere.  We  will  neglect  all 
variations  along  the  horizontal  plane.  If  we  assume  that  the  displaced  fluid  element 
is  adiabatically  compressed  by  the  surrounding  background  atmosphere,  we  may  write: 


i. 


p(z)  = P0(zo>  + 


Pq(2)  — 


= p0(zo)  + r? 


p(z)  = po(z) 


We  may  now  substitute  (1)  into  the  Euler  equation 


when  w is  the  vertical  velocity  field  of  the  fluid. 

If  we  expand  pQ(z)  in  (1)  about  zq  in  Taylor  Series,  we  obtain: 


Sp=  pD(z)  - po(zQ)  ^ — z (z  - zQ)  = - gp0(z  - zQ) 


In  first  order  (z  - zq) , we  may  hence  write, 

/ \ . (SP-1-1  t <$P  P 3p  1 92p  , \ 

Pq1^1  + — ' |£  „ p -1  i - •—  o -SL  + (z"z0) 

2 3z  Mo  2 3z  „2 

c pc  z 3z  z 

- o o o 


6P  3po 


„ I !f°  6P_  ^o  + 1 !%>  . z , 

p 3z  2 2 3z  p . 2 o 

o z p c Z o 3z  z 

O o o o 


. 3p  2 3p 

1 o g , , g o . . 

— - (z  - z ) - - — (z-z  ) 

p 3z  2 o p 3z  o 

O z C o z 

o o 


Hence, 


1 

Po  3z  z 


3p  2 

£ ^ + S_  (z  - 2 ) 

p 3z  2 o 

o c 


. , 3p  2 

1 3p  g o , g , . 

— ~ 9 - r + -hr  (z-z  ) 

p 3z  p 3z  2 o 

o c 


We  may  write 


+ oj  2 (z  ) (z-z  ) = 0 

Dt  B o o 


L 


j 
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where 


r g ^Po  2I' 

I ~ o 


l p0  a 


z c 


is  the  usual  definition  for  the  Brunt  frequency. 

We  may  now  transform  equation  (3)  into  a "particle"  equation  by  substituting 
the  trajectory  z = z(t)  into  equation  (3). 

Dw  ( z ( t ) ,t)  _ .. 

Dt 

The  final  "particle"  equation  becomes: 
d2  2 

— j (z  - zQ)  + 0)  (z  ) (z  - z ) = 0 . (4) 

, ^ o o o 

at 

This  is  just  the  equation  for  the  Brunt  oscillation  of  a small  fluid  element 
about  zq.  Normally  equation  (4)  is  derived  with  the  assumption  that  the  fluid  is 
contained  in  a weightless,  flaccid,  thermal  insulating  bag.  Such  a derivation  essen- 
tially assume  a "particle"  picture  right  from  the  beginning.  We  have  formulated  the 
problem  and  derived  equation (4)  from  a field  picture.  This  will  allow  us  to  generalize 
much  more  easily  to  the  case  when  the  background  is  no  longer  a simple  hydrostatic 
atmosphere . 

We  shall  now  consider  the  case  when  the  background  consists  of  a linearized 
Hine's  type  of  gravity  wave  propagating  in  an  isothermal  atmosphere.  A fluid  element 
in  equilibrium  with  the  background  will  no  longer  be  stationary  but  will  describe  a 
simple  harmonic  trajectory  given  by: 


z (t,z  ) = 

OO  U) 


sin  wt 


where  Uz  is  the  amplitude  of  the  gravity  wave  velocity  field  at  zq. 

Once  again  we  assume  that  if  a fluid  element  is  displaced  from  its  equilibrium, 
it  will  be  adiabatically  compressed  by  the  surrounding  background  which  now  consists 
of  both  the  hydrostatic  atmosphere  and  the  gravity  wave.  Thus,  we  may  write  in  place 


of  equation  (1) 
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PB(Z)-PB(Z0} 

p(z,t)  = pB(zQ(t))  + = PD  + 


5P„ 


B 2 
c 


(6) 


p(z,t)  = pB(zG(t) ) 


where 


PB(z,t)  = PQ(z)  + ApG(z,t) 


(7) 


p_(z,t)  = p ( z)  + Ap  (z,t) 

D O 


and  Ap  and  Ap  are  the  gravity  wave  density  and  pressure. 
G G 


Since  long  period  gravity  waves  generally  have  very  long  horizontal  wave  length, 
we  shall  again  assume  horizontal  stratification  and  write  the  Euler  equation  as: 

Dw  _ 1^  _3p  _ 


Dt  p 3z 

Substituting  equation  (6)  into  (7)  we  obtain: 

1 3P. 


(8) 


— -jp  + g = 

P dZ 


B 


6pB  3z 

PB  + 


+ g 


(9) 


1 6pB 


tPB  PB  c2  ”B 

(PB  + — ) 


6P, 


3Pr 


3z 


+ g- 


Neglecting  quadratic  terms  in  6P  , we  obtain 

B 


1 3p  ^ 

— + g 'V, 

p 3z 


1 *PB 


,_PB  3z 


+ g 


, 3P„ 

1 B B 


PB  C 


2 3z 


(10) 


1 


j 


J 


, ••  3U 

Where  ZG  = K 


G 

2 _ r i 3pb 
°h  _9Lpb  3z 


2 

c 


^ to  {1  + K X sin[oj(t  - t ) - K (z  - z ) ] } 
B z o o z G o 


(15) 


U, 


where  X = 

O 0) 


Wg  = Brunt  frequency  at  zq  . 

Equation  (14)  is  similar  to  equation(4)  except  z replaces  z and  w replaces 

yj  OH 

0)_.  Equation  914)  is  a Hills  type  of  equation  with  u>  the  gravity  wave  frequency  as 
the  fundamental  frequency  for  parametric  excitation. 

In  deriving  equation  (14),  we  have  assumed  the  Hines'  linearized  approximation. 


3 (ApG) 
3z 


+ gApG 


3U 
z 

3t 


where  Uz  is  the  Hines  gravity  wave  velocity.  According  to  Hines  (1960) , the  most 
important  non-linear  effects  arise  from  the  velocity  gradient  terms  (U  • ^)U  which 
is  neglected  in  the  above  approximation.  A more  accurate  equation  would  then  be 
given  by: 

3(AP„)  3U  DU 

— + 9ApG  = - 3T  " (U  ‘ ^)Uz  = ' -5T 
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where  AP  and  Ap  are  no  longer  the  Hines ' solution  but  each  must  now  contain 
G G 

some  correction  terms . 

Equation  (14)  would  now  be  replaced  by: 


d (z  - zG) 


+ uh2(zg)  (z~zg)  = - \ (zG)2 

c 


The  solutions  to  equation  (16)  and  other  discussions  are  given  in  the  enclosed 
preprint  (Tuan  et  al  (1976)). 

With  the  time  dependent  frequency  to  given  by  equation  (15)  , it  can  be  shown 
(Hodges  (1967))  that  so  long  as  KzXq  < 1 instability  will  not  occur  and  Richardson's 
number  will  always  be  greater  than  y . Actually,  we  are  in  fact  quite  safe  since 
Hines  (1974)  has  shown  that  even  the  usual  lower  limit  for  stability  given  by  R^  n,  i- 
is  somewhat  conservative. 

As  shown  in  Fig.  3 and  4 in  the  enclosed  preprint  (Tuan  et  al  (1976)),  "self- 
excitation" can  occur  even  in  the  presence  of  considerable  viscous  damping.  The 
energy  for  the  excitation  of  the  natural  oscillations  is  actually  fed  in  from  the 


gravity  waves. 
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FIGURE  LEGENDS 


Fig.  1.  Temperature  profiles  for  the  Summer  Season  with  exospheric  temperature 
ranging  from  600°K  to  2200°K. 

Fig.  2.  "Potentials"  for  long  period  gravity  wave  and  an  exospheric  temperature 
of  600°K.  The  dashed  curve  is  for  the  Summer  while  the  solid  curve  is 
for  Winter. 

Fig.  3.  "Potentials"  for  long  period  gravity  wave  and  an  exospheric  temperature 
of  1800°K.  The  dashed  and  solid  curves  are  for  the  Summer  and  Winter 
Seasons  respectively. 


Fig.  4. 


"Potentials"  for  long  period  gravity  wave  and  an  exospheric  temperature 


of  2200  K.  The  dashed  and  solid  curves  are  for  Summer  and  Winter  Seasons 
respectively. 

Fig.  5.  "Potential"  for  a gravity  wave  with  a 2-hour  period  and  an  exospheric 

temperature  of  1500°K.  The  figure  also  shows  the  two  resonance  states 
2 2 2 

corresponding  to  X =10.43  sec  /km  (v  , = 309.7  m/sec)  and 

phx 

2 2 2 

X = 14.33  sec  /km  (v  ^ = 264.2  m/sec).  These  two  states  are  the  only 
guided  or  partially  guided  modes  this  "potential"  will  permit. 

2 

Fig.  6.  Plot  of  phase  shift  as  a function  of  the  eigenvalue  X . The  phase  shift 

2 2 2 

goes  through  a very  sharp  change  of  tt  ratians  at  X = 10.43  sec  /km  and 
14.33  sec2/km2. 

Fig.  7.  Plot  of  the  eigenfunction  ip  for  the  Lamb  (S^)  mode. 


Plot  of  the  eigenfunction  for  the  S2  mode. 


Fig.  8. 
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Fig.  9a.  Analytic  temperature  profile  with  a very  sharp  temperature  drop  at 
the  base  of  the  thermosphere. 

Fig.  9b.  The  "potential"  for  the  temperature  profile  of  Fig.  9a  showing  a 
potential  well  with  a minimum  at  about  140  Km. 

Fig.  10.  Gravity  wave  with  a two-hour  period  and  a horizontal  phase  velocity 
of  707  m/sec. 


Fig.  11.  F-region  gravity  wave  being  ducted  by  the  earth's  curvature.  The 

strong  forward  tilt  of  the  phase  front  is  also  shown  in  the  figure. 


Fig.  12.  Comparison  of  different  gravity  wave  states  showing  the  and  S2  partially 
guided  modes  as  well  as  the  "free"  gravity  wave  modes.  (1)  is  the  F- 
region  mode  with  a horizontal  phase  velocity  of  707.1  m/sec  while  (2) 
is  a lower  altitude  "free"  mode  with  a horizontal  phase  velocity  of 
213.2  m/sec. 

Fig.  13.  "Potential"  for  a gravity  wave  with  only  a 7.5  min.  period. 
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| "ON  GRAVITY-WAVE  INDUCED  NATURAL  OSCILLATIONS" 

Abstract : 

Airglow  oscillations  (mostly  5577  8 oxygen  emission)  which  show  a 
small  period  (8-14  rnin.)  oscillation  superimposed  or.  a long  period  oscil- 
lation (several  hours)  have  frequently  been  observed.  Some  data  shows  that 
the  short  period  oscillation  is  "excited"  by  the  long  period  oscillation  as 
the  latter  increases  in  amplitude.  The  effects  often  occur  in  the  absence 
of  impulsive  energy  sources.  We  propose  a model  in  which  the  short  period 

* 

oscillation  is  relatively  localized  and  is  associated  with  the  natural 
(Brunt-Vaisala)  oscillation  of  the  atmosphere.  This  oscillation  is  trigger- 
ed off  by  a long  period  gravity  wave  through  non-linear  coupling. 

’ 


Introduction: 
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There  exists  a considerable  body  of  airglow  data  showing  a small  (6-14  min.) 
period  oscillation  superimposed  over  a long  period  oscillation  of  the  order  of  hours. 

Based  on  assuming  an  inviscid,  lossless  and  isothermal  atmosphere,  theoretical 
calculations  have  shown  that  an  impulsive  energy  source  such  as,  for  instance,  a 
bomb  explosion  can  produce  at  far  distances  (say  of  the  order  of  thousands  of  kilo- 
meters) a long  period  gravity  wave  together  with  a short-period  acoustic  frequency 
wave  (Row,  1967,  Liu  and  Yeh,  1971).  However,  the  result  requires  the  short  period 
acoustic  frequency  wave  to  arrive  first  followed  by  the  long  period  wave. 

Much  of  the  airglow  data  mentioned  in  the  first  paragraph  are  not  related  to 
any  known  impulsive  energy  sources.  Furthermore,  some  of  the  data  appears  to  shew 
that  the  short-period  oscillation  actually  follows  rather  than  procedes  the  arrival 
of  the  long-period  oscillation  and  is  slowly  excited  as  the  long  period  oscillation 
increases  in  amplitude. 

t We  propose  to  examine  a model  (Tuan  et  al.,  1975)  in  which  we  consider  the 

short-period  oscillation  to  be  local  and  to  be  associated  with  the  natural  Brunt- 
Vaisala  oscillation  of  the  atmosphere.  The  long  period  oscillation  will  be  assumed 

i 

to  ce  a gravity  wave  which  may  bo  generated  a long  distance  away.  The  non-linear 
part  of  the  gravity  wave  serves  as  a forcing  term  which  couples  the  B-V  oscillation 
with  the  gravity  wave.  In  the  limit  of  small  oscillation  amplitudes,  the  gravity 
wave  becomes  linear  and  we  have  essentially  damped  simple  harmonic  motion.  As  the 
amplitude  of  the  gravity  wave  increases,  departure  from  linear  behaviour  occurs. 

If  the  amplitude  of  the  gravity  wave  increases  just  slightly  beyond  a certain  value, 
the  non-linear  term  becomes  strong  enough  to  excite  natural  oscillation  and  we  have 
a superposition  of  both  the  natural  oscillation  and  the  gravity  wave  oscillation. 


I 


Summary  of  Observed  Results: 

A typical  superposition  of  a small  period  (s,  8-10  min.)  oscillation  superim- 
posed over  a long  period  oscillation  for  the  5577  m 01  is  shown  in  Fig.  1 (Okuda, 
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1975) . The  data  is  taken  at  Sendai  in  January  1966. 

Fig.  2(a)  (Okuda,  (1962),  Silverman  (1962))  shows  the  interesting  case  where 
the  long  period  oscillation  increases  in  intensity  producing  the  "self  excitation" 
in  the  short  period  oscillation.  To  a much  lesser  extent,  the  Oct.  28-29  (1961) 
result  of  Fig.  2(b)  also  shows  this  effect.  By  contrast,  on  a relatively  calm 
day  (Oct.  21-22)  there  is  relatively  much  less  oscillation.  The  very  gradual  nature 
of  the  introduction  of  the  short  period  oscillation  which  essentially  follows  rather 
than  precedes  the  long  period  one  seems  to  suggest  that  the  excitation  is  local  and 
most  certainly  not  impulsive. 


Theory : 

In  the  absence  of  gravity  waves,  a vertical  displacement  of  a fluid  element  will 
immediately  bring  about  a restoring  buoyancy  force.  For  small  displacements,  neglec- 
ting viscosity,  simple  harmonic  motion  is  produced.  It  is  given  by: 

d2  U - z ) 

— — f + vzo)(z-zo)  = 0 (1) 


where 


-g 


tY-ng2  + 

2 

c 


Po  = pq(z)  is  the  background  density 

zq  = the  position  at  which  the  fluid  element  is  in  equilibrium  with 
the  background. 

Equation  (1)  may  be  obtained  directly  from  the  Euler  equation  if  we  assume 

Vz)  “Po(zo) 

p (z)  = p (z  ) + (2) 

o o 2 

c 

P(z)  = P (z)  = background  pressure  (3) 

Equation  (2)  gives  the  density  of  the  fluid  element  when  it  is  displaced  slightly 
from  zQ  to  z and  adiabatically  compressed  lay  the  background  pressure  change. 

In  the  presence  of  a gravity  wave,  we  may  assume  that  the  background  now  con- 


sists of  both  an  isothermal  hydrostatic  atmosphere  and  the  gravity  wave. 


'hus , 
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we  may  write 

PB(z,t)  = PQ(z)  + APG(z,t)  (4) 

PB(z,t)  = PQ(z)  + ApG(z,t)  (5) 

where  AP  and  Ap_  are  the  gravity  wave  pressure  and  density.  If  a fluid  element 
G G 


is  displaced  from  zG(t)  to  z,  its  density  becomes: 


p (z,t)  = pB(zG(t) ) + 


P(z,t)  = PB(zG(t)) 


Pb(z)-Pb(zg) 


where  z£  = z (t,z  ) is  the  trajectory  of  the  fluid  element  if  it  remains  in  equili- 


brium with  the  background.  Substituting  equation  (4) , (5) , (6)  and  (7)  into  the  Euler 
equation,  we  obtain,  after  suitable  expansion,  the  following: 
d2(z-z  ) 

■= + u>  (z  ) (z  - z ) = 0 (8) 

dt2  H G G 


where 


2 . . 
UH  zg^ 


i apT 


pb  3z 


'U  w {1  + k X sin [w (t  - t ) - k (z  - z ) ] } (9) 

B Z O O Z G O 

In  equation  (9) , k^  is  the  vertical  wave  number,  zG  = Xq  sin  ait  and  Xq  is  the  ampli- 

2 

tude  of  the  gravity  wave  oscillation.  It  is  easy  to  show  that  wH  given  by  equation 

2kzUT0 

(9)  is  identical  to  the  expression  given  by  Hodges  (1967)  for  — >>  1.  We  will 

consider  for  the  rest  of  this  paper  the  special  case  where  kzXQ  < 1-  This  means 
that  the  Richardson's  number  will  always  be  greater  than  1/2  so  that  instability 
will  not  occur,  (Hodges,  1967). 

In  deriving  equation  (8)  we  have  assumed  the  Hines'  linearized  approximation. 


3(APC) 


+ gApG  = " 3t~ 


where  Uz  is  the  Hines  gravity  wave  velocity. 


According  to  Hines  (1960) , the  most  important  non-linear  effects  arise  from 
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the  velocity  gradient  term  (U  • V^)U  which  is  neglected  in  the  above  approximation. 
A more  accurate  equation  would  thus  be  given  by: 


8(APg) 


+ g ApG  = - TT 


(U.V)U2  = ___ 


where  AP  and  Ap  are  no  longer  the  Hines ' solution  but  each  must  now  contain 
G G 

some  correction  terms. 

Equation  (8)  would  now  be  replaced  by: 
d^  (z  — z ) 

- ~ + u 2 (z  ) (z  - z ) = (U‘V)U=-~  (z  ) 2 (10) 

2 Ho  O Z ZO 

dt  c 

Equation  (10)  reduces  to  equation  (8)  for  small  gravity  wave  amplitudes.  We  would 
like  to  emphasize  that  while  we  have  not  specifically  included  viscosity  effects 
in  our  discussion,  we  did  include  them  phenomenologically  in  our  calculations. 

They  are  not  important  and  serve  merely  to  damp  out  natural  oscillations  arising 
from  transients.  We  shall  leave  them  together  with  a more  complete  derivation  of 
equation  (10)  for  later. 

Results : 


Figures  (3)  and  (4)  are  the  theoretical  results.  We  have  assumed  a Brunt 
Period  of  4.9  minutes  and  an  isothermal  background  atmosphere.  Figure  (3)  shows 
what  happens  after  the  transients  have  been  damped  out  by  viscosity.  The  amplitud 


XQ  of  the  gravity  wave  is  1.58  Km.  Even  with  this  much  amplitude,  the  Richardson's 
number  is  still  greater  than  1/2  so  instability  will  not  occur.  However,  some 
significant  "self-excitation"  of  the  natural  oscillation  due  to  the  non-linear  fore 
ing  term  occur.  The  period  of  such  oscillations  is  of  the  order  of  5 minutes 
close  enough  to  the  B-V  oscillation. 

Figure  (4)  shows  what  happens  when  we  use  heavy  viscous  damping.  In  general 
viscous  damping  depends  on  the  size  of  the  displaced  fluid  element.  The  smaller 
the  fluid  element,  the  greater  the  damping  effect  becomes.  In  Fig.  (4)  we  see  tha^ 
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the  transient  effects  are  damped  out  in  about  one  cycle  and  the  "self  excitation" 
begins  after  about  3 cycles  of  the  gravity  wave. 

Discussion 

While  we  do  not  claim  that  the  superposition  of  a small  period  oscillation  of 
the  order  of  B-V  frequency  and  a longer  period  gravity  wave  oscillation  is  always 
caused  by  non-linear  coupling,  the  general  shape  of  the  oscillation  curves  and  the 
manner  in  which  the  natural  oscillation  is  excited  do  suggest  that  this  may  well 
be  a plausible  mechanism.  We  wish  to  emphasize  that  it  is  not  surprising  that 
"self  excitation"  can  occur  since  equation  (8)  which  is  rigorously  derived  is 
essentially  a Hills  type  of  equation  whose  solution  can  be  unstable.  In  fact,  any 
disturbance  introduced  into  equation  (8)  would  produce  natural  oscillations  which 
amplify  with  time  and  is  only  prevented  from  rapid  increase  by  viscous  damping. 

The  energy  is  fed  into  the  natural  oscillation  from  the  gravity  wave.  Unlike 
equation  (8) , the  derivation  of  equation  (10)  involves  a number  of  assumptions  and 
must  hence  be  considered  as  semi-phenomenological.  All  we  wish  to  show  here  is 
that  a (U  • V)U  type  of  non-linearity  is  sufficient  to  generate  natural  oscilla- 
tions  even  when  Richardson's  number  is  still  relatively  large. 


\ 
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Figure  Captions 


Figure  1:  The  5577  8 01  emission  is  plotted  against  local  time.  The  data 

is  taken  at  Sendai  by  Okuda. 

Figure  2:  (a)  This  plot  of  5577  8 01  emission  provides  a good  example  showing 

how  a small  period  oscillation  may  be  slowly  excited  by  the  long 
period  gravity  wave.  (b)  This  shows  a similar  but  less  obvious 
self  excitation  at  Sacremento  Peak.  A comparison  is  also  made  with 
a quiet  day. 

Figure  3:  A small  amplitude  gravity  wave  with  a one  hour  period  is  shown  oscil- 

lating about  95  km.  The  wave  approaches  a sinusoidal  shape  for  small 
enough  amplitudes.  The  natural  oscillation  is  damped  out  since  for 
small  amplitudes  equation  (10)  reduces  to  equation  (8) . 

Figure  4:  Gravity  wave  with  a period  of  one  hour  but  with  10  times  the  damping 

(a  much  smaller  fluid  element)  is  shown  in  Fig.  4.  The  transient 
natural  oscillation  is  damped  out  after  one  hour.  Self  excitation 
of  natural  oscillation  began  after  3 hours  with  much  less  oscillation 
amplitude  than  shown  in  Fig.  3.  The  present  example  is  an  extreme 
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While  the  experimental  observations  of  airglow  intensity  fluctuations 
end  iis  possible  connections  with  TJD  have  been  made  by  i:.any  investigators, 
very  little  attention  has  been  given  to  a theoretical  analysis  of  such  ion- 
ospneric  disturbances  on  airglow.  In  this  paper  wc  specifically  examine  the 
effect  of  an  atmospheric  gravity  wave  on  the  6300  ?>  01  end  the  5200  % KI 
emission  in  mid-latitudes.  Yte  will  also  shov.-  that  for  the  6300  ?.  01,  if  we 
assume  reasonable  values  for  the  ionospheric  para.  :terr.  as  veil  as  a gravity 
wave  with  a period  of  approximately  3 hours,  the  resultant  fluctuations  in 
airglow  is  consistent  with  Dach's  observed  data  as  well  as  with  the  simultan- 
eously measured  variations  in  the  ionospheric  parameters. 


Introduction 

An  association  of  airglow  intensity  variations  with  traveling  ionospheric 
disturbances  has  been  established  by  Barbier  (1964,  1965)  for  the  6300  a 01 
emission,  and  by  Weill  and  Christophe-Glaume  (1957)  for  both  the  6300  $ Cl  and 
5200  n NI  emissions  by  a direct  comparison  of  the  two  types  of  measurements. 

In  addition,  some  published  data  (Silverman,  1962;  Okuda,  1362;  Dachs,  1353) 
can,  in  retrospect,  be  seen  as  having  been  probably  due  to  similar  causes. 

The  data  of  Dachs  include  observations  on  both  5577  ?\  01  and  6300  01  as  well 

as  of  the  ionospheric  parameters  foF2  and  h'F2,  and  are  thus  particularly 
valuable  for  analysis  from  the  point  of  view  of  a possible  wave  origin  of  the 
variations.  Disturbances  in  airglow  can,  of  course,  arise  from  a variety  of 
causes,  and  simultaneous  measurements  of  other  parameters  are  needed  to  establish 
the  cause  of  a disturbance. 

Over  the  past  few  year’s  considerable  theoretical  and  experimental  work 
has  been  carried  out  on  the  effects  of  atmospheric  gravity  waves  on  ionospheric 
parameters.  The  effects  of  these  waves  on  optical  emissions,  however,  has 
received  no  attention  from  a theoretical  viewpoint.  A suitable  analytic  model 
for  such  a treatment  has  been  presented  recently  by  Porter  and  Tuan  (1973). 

This  model  for  the  response  of  the  F-layer  to  gravity  waves  differs  fundamentally 
from  the  usual  theoretical  treatment  in  several  respects.  Of  these  the  most 
pertinent  for  airglow  applications  are:  (1)  The  analytic  model  is  fully  time- 

d ?3> 3nd.tiTi^‘  • Cr.G  IT 3 V 2.  +1  Q ^ 3 LI C 1*1  3.  T?.OCl.0.1  1.S  rLV*?n  ?.  3 ir.USOi^  3^.  cr^av*  ! ** 


wave,  the  response  of  the  electron  density  oscillation  is  no  longer  in  general 


sinusoidal  and  the  time  behavior  varies  with  height.  This  makes  it  more  suitable 


B-4 


for  studying  both  the  general  decay  of  the  F-layer  and  airglow  at  night  as  well 
as  such  detailed  behavior  as  differences  in  phase  between  the  oscillations  of 
electron  density  and  the  airglow.  (2)  The  analytic  model  remains  valid  even  if 
the  percentage  oscillations  in  the  electron  density  or  airglow  is  large  compare 
with  ambient  values.  This  is  especially  important  for  the  6300  ° 01  emission 
whose  luminosity  profile  peaks  at  a position  (about  one  scale  height  below  F2 
peak)  where  the  electron  density  and  airglow  oscillations  are  rather  large. 

In  this  paper  we  first  outline  the  computation  for  the  6300  S 01  emission 
intensity  during  the  passage  of  a gravity  wave.  For  long  period  gravity  waves 
the  continuity  equation  from  which  the  emission  intensity  of  the  6300  a 01  can 
be  calculated  reduces  to  the  quasi  equilibrium  theories  of  Peterson,  Van  Zandt 
and  Norton  (1066),  in  which  the  production  of  the  excited  0(~D)  states  is  as sum 


d 


2'd 


to  equal  the  loss  through  optical  emission  and  quenching.  Comparison  of  the 
theoretically  calculated  changes  in  both  electron  density  and  luminosity  pro- 
files will  serve  to  bring  out  the  essential  features  of  the  theory. 

Vie  shall  next  consider  the  time  behavior  of  the  total  electron  content  as 


compared  with  the  airglow  oscillations.  It  will  be  shown  that  the  total  electron 
content  shews  less  oscillatory  behavior,  and  that  only  by  plotting  the  time 
derivative  of  the  total  electron  content  can  the  oscillations  be  more  significant 
However,  such  fluctuations  are  detectable  (Davis  and  daRosa,  1969).  The  reduced 
oscillation  in  the  total  electron  content,  coupled  with  vine  large  otical  emission 
oscillations,  provides  us  with  a possible  criteria  for  distinguishing  between 


Long  pari.oi  gravitj  nave  effects  and  those  pr; 


:ea  _-y  suen  ocr.er  sources  suc.o 


as  particle  participation. 
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We  will  then  consider*  the  effect  of  a gravity  wave  on  the  5200  9 NI  emission 


intensity.  Here  we  can  no  longer  use  the  equilibrium  theory  that  we  have  been 

using  for  the  6300  9 line,  since,  at  least  at  higher  altitudes,  the  production 
2 

of  N(  D)  states  can  no  longer  be  balanced  by  either  the  radiative  loss  or  the 
quenching.  We  will  also  consider  the  importance  of  electron  quenching,  which  is 
not  significant  for  the  6390  9 01  line,  but,  as  we  shall  see,  is  much  more 
important  for  the  nitrogen  line. 

The  theory  derived  here  will  then  be  used  for  the  interpretation  of  the  experi- 
mental results  of  Dachs  (1963)  at  Tsuraeb,  South  West  Africa.  These  include  simul- 
taneous measurements  of  foF2,  h'F2,  and  the  airglow  intensity  variations  with  time, 
thus  allowing  for  a test  of  both  aspects  of  the  theory. 

Finally,  we  shall  also  study  the  effect  of  varying  the  quenching  rates  on  the 
phase  of  the  oscillations.  We  will  consider  this  phase  difference  in  relation  to 
the  effect  of  other  mechanism  which  can  also  cause  a phase  difference. 
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Tims  behavior  of  electron  density  and  6300  ?!  01  luminosity  profiles 

Theoretical  treatments  of  the  630C  n 01  emission  taking  account  of  different 
factors  and  calculating  various  parameters  are  available  from  a number  of 
authors  (see,  for  example,  Chamberlain,  195C;  Lagos,  belle;-;  and  Silverman,  1S63: 
Lagos,  1954;  Peterson,  Van  Zandt  and  Horton,  1966;  and  Tuan,  1959).  The 
following  treatment  leans  most  heavily  on  the  work  of  Paterson,  Van  Zandt  and 
Norton. 

We  make  the  following  assumptions:  (for  their  justification  see  Peterson, 

Van  Zandt  and  Norton  1366);  (1)  the  time  derivative  term;  in  the  continuity 
equation  for  the  density  of  oxygen  atoms  in  their  Q(^D)  states  may  be  neglected; 

(2)  the  diffusion  term  is  neglected;  and  (2)  we  assume  chemical  equilibrium  for  the 

• "f"  -i-  _j 

production  and  loss  of  C2  molecules  as  well  as  neglecting  the  0^  and  NO  molecule 
in  the  equation  for  charge  neutrality.  No  then  obtain  the  following  well  known 
expression  for  the  volume  emission  rate  c , for  the  6300  2 lire, 

- f.  ir.n  5 


e6300  (2’t)  = n 


IA620o)  y [,  02] 
l^D  j 1 + 4— 


<S1  CK2]  + S2Kc> 


(1) 


where  >/  is  the  efficiency  for  recombination,  y is  the  rate  coefficient  for  charge 

*f* 

transfer  between  0 and  0?  , S and  S2  are  the  nitrogen  and  electron  quenching 
coefficients  respectively.  0oh  profiles  obtained  by  Sckunk  and  Walker  (preprint) 
could  later  be  used  in  equation  (1).  However,  since  the  same  chemical  equilibrium 

and  charge  neutrality  ipproximat ion s have  beets  used  in  the  xpressicn  rcr  f , 


(Porter  and  Tuan,  1974),  we  have  decided  to  maintain  consistency  and  use  equation 
(l)  for  the  emission  rate.  Since  our  primary  interest  is  in  the  basic 
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dynamic  response  of  the  airglow  to  gravity  waves  rather  than  any  detailed 
behavior,  ue  have  also  neglected  the  contribution  from  dissociative  recombination 
of  N0+  (Silverman,  1970)  as  well  as;  cascading  from  the  "^S  states  in  the  production 
of  D states  (Rishbeth  and  Carriott,  1969).  V.Te  will  also  neglect  the  time  varia- 
tion in  0^  and  If,  caused  by  the  gravity  'waves.  In  reality  the  equations  for  such 
neutral  atmospheric  constituents  are  really  coupled  to  the  equation  for  the  ion- 
electron  plasma.  To  simplify  all  these  processes  we  have  used  stationary  models 
for  0^  and  In  practice,  the  approximations  should  be  good  since,  as  already 

mentioned,  the  percentage  variation  in  the  neutral  gas  caused  by  the  gravity  wave 
is  usually  very  much  less  than  the  percentage  variation  for  the  ion -electron 
plasma,  especially  in  important  regions  such  as  the  peak  of  the  luminosity  profile. 
By  fur  the  most  important  contribution  to  the  oscillation  in  the  airglow  arises 

from  changes  in  the  electron  density  Ke  in  equation  (1). 

1 

Host  recent  theoretical  treatments  for  the  time-dependent  behavior  of 
F- layer  elect!  >n  c nsity  (e.g.  Testud  and  Francois,  1971;  Klostermeyer,  1972  a,b) 
under  the  irfluance  of  gravity  vraves  are  essentially  based  on  perturbation 
theory  in  vjhich  the  unperturbed  ionosphere  is  assumed  to  be  a stationary  Chapman 
function  independent  of  time.  Also,  in  perturbation  treatments  the  tin.?  depen- 
dence for  the  p ’ is  assumed  to  be  sinu  oidal  if  the  gravity  wave  is  sinusoidal. 


and  a simple  product  form  for  the  time  and  spatial  dependence  is  usually  assumed. 

The  luminosity  peak  for  the  night-time  6300  ?.  line  occurs  about  one  scale 
height  below  the  F2  peak  (Tuan,  1959)  where  the  electron  density  oscillation  : -r. 
be  larger,  (Porter  and  Tuan,  1974).  As  an  ux ample , for  the  gravity  wave  that  Thome 
(1968)  analysed,  the  oscillation  in  electron  density  at  this  height  would  be  between 
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1 


where  z ='  h-hQ  is  the  difference  between  the  actual  height  h and  some  arbitrary 

reference  height  h ; !!  is  the  scale  height;  p is  the  scale  height  factor; 

Ao=Ao(p)  is  a decay  constant  (see  Appendix);  w is  the  angular  frequency  of  the 

gravity  wave;  q is  a parameter  (0  <_  q <1/2 ) which  controls  the  damping  of  the 

gravity  wave  (when  q = 1/2  we  obtain  the  undamped  Hines  Model  (Hines,  I960); 

B(q)  is  a coefficient  depending  on  gravity  wave  and  ionospheric  parameters 

(see  Appendix);  t is  some  initial  time  giving  us  the  phase  of  the  gravity  wave 
o , — - 


at  t = tQ;  a = 2H 


Pp 


_ where  6 and  D are  the  recombination  ana  diffusion 
D o o 

o 


coefficients  at  hQ;  0 is  the  angle  between  the  magnetic  field  and  the  zenith; 

A is  the  amplitude  of  the  gravity  wave  at  hQ;  ln(p,w,t)  and  R (p,u,t)  are 

oscillating  functions  of  time  defined  in  the  appendix  F (p)  and  G (q)  arc 

factors  which  control  the  dependence  on  the  scale-height  factor  end  the  damping 

. -3 /2  . 

of  the  gravity  wave  respectively  and  are  given  m the  appendix;  L ' is  a 
Laguerre  Polynomial  of  half  integral  order. 

We  see  at  once  from  equation  (2)  that  if  we  ignore  the  contribution  from 
the  infinite  sum  on  the  right  hand  side,  we  would  have  a Chapman  function 
multiplied  by  an  exponential  damping  together  with  oscillations.  This  is  as 
would  be  expected  since  for  the  night-time  situation  considered  here  production 
can  be  excluded.  (We  may  easily  include  production  if  required). 

The  infinite  sum  on  the  right  hand  side  of  equation  (2)  reduces  to  a few 
terms  for  a gas  either  perfectly  mixed  or  in  diffusive  equilibrium  and  for  a 
constant  velocity  amplitude  (independent  of  height)  for  the  gravity  wave  (q  = 0), 
(see  Appendix). 


position  (Fig.  1(c)).  If  we  had  plotted  the  luminosity  profiles  without 
quenching,  the  peak  luminosity  would  have  moved  a considerable  distance  (about 
to  Ka)  downward  when  the  12  peak  comes  down.  This  implies  that  xhe  luminosity 
peak  is  essentially  kept  up  by  the  rapid  increase  in  quenching  at  lower  altitudes 
The  reason  that  it  would  only  move  15  km  rather  than  the  20  km  for  the  F2  peak 
is  due  to  the  rapid  drop  in  available  0^'  for  the  production  of  states  at 
higher  altitudes  which  tends  to  keep  the  luminosity  peak  down.  Thus,  the  com- 
bined effect  of  both  is  to  keep  the  luminosity  peak  from  having  larger  oscilla- 
tions in  xhe  vertical  direction. 


Oscillation 


otal  electron  content  and  6300  A 01 


_ .entity  for  the  G20C  h 01  v:c  integral 

equation  (1)  ever  height.  The  results  are  plotted  against  local  tine  in  Tig. 
We  have  also  plotted  the  total  columnar  electron  content  as  well  as  the  nag.-: 
time  derivative  of  this  total  electron  content.  The  parameters  used  here  a 
the  same  as  those  adopted  for  Fig.  1.  One  general  feature  is  that  all  throe 
curves  decay  with  time,  an  expected  result  since  no  ion-pairs  production  r.ech 
have  been  included.  The  most  significant  feature  is  that  while  there  is  con- 
siderable oscillatic wi  a time  in  the  airglow  intensity,  the  total  electron 
content  shows  relatively  little  fluctuation.  The  lack  of  fluctuation  in  -he 
total  electron  content  is  readily  unde? standable . We  are  first  of  all  dealin 
with  a night-time  situation  when  v:e  may  neglect  photoprcduction.  The  flue tun 
tions  arise  when  the  ions  are  moved  alt<  rnatoly  into  rarified  ; tme sphere  wher 
the  recombination  is  less  and  denser  atmosphere  where  the  recombin  tion  is  gr 
Thu  fluctuation  is  further  inhibited  by  phase  cancellation  between  the  i oti  n 
lens  above  and  below  the  F2  peak.  Since  we  a <’  aling  with  a horizontally 
ducted  neutral  wave  with  no  vertical  pi....:-  pro]  ress  ion,  the  phase  variation  . 
ions  depends  almost  entirely  on  the  shaj  ■ of  the  F-lay>  r.  Fig.  6 shows  how  . 
ph  e of  the  ion  motion  above  the  vicinity  of  the  F2  peak  is  in  general  opj 
to  that  below  th(  peak.  We  have  I 23  2 an  example  of  po:  siblc  ph  s<  cancel!  . r 
even  when  the  gravity  wave  docs  not  pc'uoss  phase  variation  along  the  line  of 
' ' * C ' \'r  V 1 ( l • >')  " • y . . ’ . ‘ . * '•j  . r.~  \"L  . * . r.  ' **  z, ; iMT1 


from  phase  variations  in  the  gravity  wave  itself. 


B-13 


The  columnar  airglow  intensity,  on  the  other  hand,  is  proportional  to  the 
ion-pair  recombination  if  we  neglect  quenching.  By  integrating  the  continuity 
equation  for  electrons  over  z and  dropping  the  production  tern  we  obtain  the 
result  that  the  columnar  recombination  rate  is  equal  to  the  negative  of  the  rate 
of  change  of  columnar  electron  content  since  the  divergence  term  in  the  continuity 
equation  drops  out  from  Gauss  Theorem.  Thus,  if  the  quenching  v.-erc  neglected, 
we  v.’ould  expect  the  airglow  intensity  to  he  directly  proportional  to  the  negative 
of  the  rate  of  change  of  total  electron  density.  This  is  shown  very  clearly  in 
l’ig.  2 in  which  the  airglow  fluctuations  follow  the  negative  of  the  rate  of  change 
of  total  electron  content  closely.  Both  curves  are  opposite  in  phase  to  the  total 
electron  content  as  expected.  If  the  quenching  terms  are  very  large,  then  a phase 
difference  will  be  introduced  between  the  airglow  fluctuations  and  the  rate  of 
change  of  total  electron  content.  In  general,  there  are  two  possible  ways  to 
introduce  this  phase  difference.  One  is  to  introduce  quenching  as  we  have  just 
mentioned.  The  other  is  to  introduce  some  production  process.  This  point  will 


be  further  discussed  later. 
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The  5200  n Nitrogen  Line 


The  5200  HI  is  produced  in  the  F-regicr.  by  dissociative  recombination  in  a 


LiiJii'jcr  r j-!..-. _ i’.yy  7.^ 


f rccj 


u il  o-  tro  c ■ : cccs  ci  c_  r,':  vv  vc\  c c 


this  line  is  cf  special  interest  because  of  the  very  long  lifetime  of  the  nitrogai 


“Dg/g^tates . We  con  r.c  longer  use  the  quasi-cquilibri u;r.  theories,  especially 


higher  altitudes  where  quenching  can  no  longer  compensate  for  the  very  slew 


radiative  less.  If  we.  neglect  diffusion  (Hernandez  and  Turtle,  1262),  the 


2 j 

continuity  equation  for  N(  D),  the  atomic  nitrogen  density  in  the  D state  is 


given  by: 


?-N(  2D) 


= p - (a2„  + K1 + 


In  equation  (3)  P is  the  production  rate  for  H(  'D),  and  K„  are  the  quenching 


coefficients  for  0-  and  electrons  respectively,  A„  is  the  Einstein  coefficient  for 

2 * D 
h(  D).  If  we  go  down  to  low  enough  altitudes  where  the  quenching  is  sufficiently 


strong  to  balance  the  production  rate,  we  may  neglect  the  time  derivative  term  and 


the  theory  reduces  to  a quasi-equilibrium  theory  exactly  analogous  to  that  of 


Feterson,  et.al.  (1266).  For  the  production  r«.te  v.e  again  assume  that  the  pro- 

duct ion  cf  1!(  D)  through  dissociative  recombination  of  HO  is  balanced  by  charge 

+ 2 
transfer  between  0 and  N . In  this  way  the  production  rate  of  N(“D)  states 


varies  with  height  as  the  product  of  neutral  11^  concentration  and  H . 


The  loss  rate  of  K(  D)  is  represented  by  the  three  terms  on  the  right-hand 


side  of  equation  (3),  corresponding  to  radiative  loss,  0o  quenching  and  quenching 


due  to  superelastic  collisions  with  the  electrons.  We  have  assumed  that  A = 1/26  hr 
(Chamberlain,  1261)  and  = 8 x 10  ^ cm^  sec  ^ (Seaton,  1256).  For  the  rate 
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coefficient  Ik  with  respect  to  quenching  v;e  have  tried  both  = 1.4  x 10 
cm°  sec  (Slanger,  et.  al. , 1971)  and  = 1.4  x 10  ^ cm3  sec  ~ (Hernandez 


and  Turtle , IScO).  The  two  values  differ 


order  of  ::.c 


UC  . 


is  reason  to  believe  that  the  value  given  by  H 
small  sice  (Hernandez,  private  ecri-municst  ion ) . 
integrated  subject  to  the  initial  boundary  con 


iernandez  and  Turtle  may  be  on 

Fquaticrj  (3)  may  be  simply 
2 

.ditions.  i'rorn  che  '■(.  'D)  ccn- 


the 


centration  we  may  readily  calculate  the  columnar 

4 2 

the  ( S„ - D_  ..,,.)  forbidden  optical  transi- 
a/ 2 b/2 , o/2 

The  results  arc;  sho;m  by  the  solid  curve  in  Fig. 


zenith  ai"'g-ov>'  intensity  for 
ion  as  a function  of  time. 

3 for  the  quenching  coefficient 


of  Slanger,  et.  al.  (1971)  ar.d  the  dotted  curve  for  the  coefi iciert  used  by 
Hernandez  ar.d  Turtle  (96S).  For  both  curves  re  have  used  a scale-height  factor 


of  p = 1.7b.  All  otbei  gr  ity  i iv<  parameters  are  again  the  same  an  those  us<  1 
to  fit  Theme's  results,  (Porter  and  Tuan,  loc,  cit.). 

The  three  Important  features  are  (l)  the  oscillations  for  the  higher  0^ 
quenching  rate  of  Slanger,  ct.  al.  (1971)  is  less  pronounced  than  those  pro- 
duced by  the  lower  quenching  coefficient  of  Hernandez  and  Turtle  (1909),  (2) 
there  is  a phase  advance  with  the  higher  quenching  coefficient  of  between  10 
to  20  minutes.  The  reasons  for  this  will  be  given  in  a later  section,  (3) 
the  luminosity  peak  for  the  higher  quenching  coefficient  is  in  general  about 
50  1cm  higher  than  the  peak  for  the  lower  coefficient. 

In  comparing  those  results  with  the  behavior  of  the  6300  a line,  there  are 
the  following  k.  pert ant  differences : The  electron  quenching  fer  the  nitrogen 

line  is  relatively  speaking  much  more  important  than  for  the  oxygen  line.  If 
we  use  the  coefficient  of  danger,  et.al.  (1971),  the  quenching  due  to  molecules 
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for  the  nitrogen  line  is  down  by  more  than  an  cider  of  magnitude  from  the 
quenching  for  oxygen  due  to  the  combined  effect  of  a smaller  quenching  coefficient 

t-  well  ..  smaller  abundance  of  o.  molecule  in  comparison  with  b.  molecules  . 

Z ~ -a 

The  difference  in  electron  quenching  in  only  about  a factor  of  2. 

The  relatively  more  important  electron  quenching  has  the  effect  of  holding 
down  the  £200  a luminosity  peak,  since  the  available  electrons  for  quenching 
increase  until  the  F2  peak,  well  above  the  luminosity  perks  of  both  the  COCO  8 
and  the  £200  8 lines. 

The  net  result  is.  that  although  one  ir.ay  at  first  expect  the  5200  8 iuviiuv fty 
pea]  to  I well  < ove  tl  oxygen  peak,  owing  to  a much  longei  half-li  ",  the 
r<  duced  mo]  qu  Lng  at  lovfer  altitudes  coupled  with  a r<  L<  . ivt  in< 

in  electron  ] ing  ( unimportant  for  ibe  6300  8 line)  at  higher  altitudes 

would  larii  down  the  position  for  the  5200  8 line  to  comparable  height 

levels  £ ■.  s that  for  the  . On  line.  This  is  not  in:  n ist<  t witl  the  obs<  r- 
vations  of  brill  and  Christophe-Glaume  (1S67)  who  argued  fie:,  gee:  ■ :.c  con- 

siderations that  the  5200  8 peak  is  about  12  km  above  .the  6300  8 p-  ck  if  the 
speed  of  the  assumed  T1D  is  assumed  to  be  the  same  at  the  peaks  of  both 
emissions.  Vfe  may  also  at  the  same  time  conclude  that  since  both  peaks  occur 
at  nearly  the  same  height,  the  oscillation  amplitudes  should  be  cer.p  arable . 

This  is  borne  cut  by  comparing  Fig.  3 with  Fig.  2. 
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Comparison  with  simultaneous  observations 


fcF2,  h'F2  and  6300  8 01  emission 


We  have  mentioned  earlier  that  some  published  data  seems  likely  to  have 
e:diibited  the  effects  of  gravity  waves.  In  this  section  we  use  some  cf  Bach' 
(1968)  simultaneous  measurements  of  ionospheric  parameters  and  the  6300  ? 01 
emission  to  illustrate  the  theory  by -fitting  the  observations  to  an  assumed 


gravity  wave. 

Bachs  (196S)  has  made  intensity  measurements  of  the  $300  ?%  and  the  5577  8 
emission  lines  as  a function  of  time.  These  measurements  together  with  sim.ul- 
taneous  measurements  of  fcF2  ar.d  h'F2  seem  to  show  cons  Lderabli  orrelation, 
especially  between  the  63C0  ° line  and  the  ionospheric  par am .tors.  In  Fig.  4 
vie  have  shown  a comparison  between  our  theoretical  results  and  the  experimental 
values  for  foF2,  h'F2  and  the  6300  ?,  airglow.  For  the  theoj  . Lc  f | v ■ t<  .. 
we  have  chosen  p = 1.5,  H = 65  Km,  the  recombination  coefficients  6(200  Km)  = 
0.34  hr  '"*■  and  the  diffusion  coefficient  D(200  Km)  --  6.1  >:  103  Km 3 hr  \ These 
are.  somewhat  different  from  the  parameters  chosen  to  fit  Thome's  observation  but 
considering  the  fact  that  v:e  have  completely  ignored  production  processes  (the 
observations  are  made  between  1800  and  2400  L.T.)  they  are  the  right  order  of 
magnitude.  Tor  the  gravity  wave,  we  assume  the  velocity  amplitude  A(200  Km)  = 
35.2  n see  ^ and  the  damping  factor  q = 0.49  which  is  rather  close  to  the  Hines 
(1950)  undamped  model.  The  quenching  at  200  Km  is  assumed  to  be  0.24  sec  1 and 
the  electron  quenching  is  neglected.  The  quantum  efficiency  for  mission  from 
the  ■'‘D  is  assumed  to  be  0.5. 


In  view  of  the  rather  good  agreement  between  theory  and  experiment  for  both 
the  airglow  and  the  ionospheric  parameters,  one  cun  readily  conclude  that  even 
without  a total  electron  content  measurement  we  have  a gravity  wave  passing 


A 
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through  with  a period  of  between  2 ■ 
airglow  intensity  and  the  h'F2  show 
v;hc-n  the  airglow  intensity  is  at  a r 


to  3 hours.  The  comparison  between 
all  the  expected  phase  correlations 
r..inar.;um  vr. a i'2  peak  is  upproxirsat ely 


the 
, i.e. 

cLu  ci 


maximum. 
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Effect  of  quenching  and  the  criteria  for  different  disturbances 

As  already  mentioned  in  a previous  section,  quenching  can  cause  a phase 


c.Titrerce  m tr.e  os 


Llation  of  ttv.  air  glow  intensity. 


K 01  for  two  Ng  quenching  coefficients.  The  solid  curve  uses  a quenching 
- . . -1]  3 -i 

coefficient  = 1.0  x 10  ' cm  sec  ' while  the  cashed  curve  uses  2^  = 2.7  x 

_J1  3 _i 

10  ' cm  sec  . All  other  parameters  are  the  same  as  those  used  to  fit  Daeh's 
(1969)  experimental  observations.  The  higher  quenching  coefficient  seems  to 
produce  less  oscillations  as  well  as  a phase  advance  by  about  half  an  hour. 

It  is  easy  to  understand  why  the  oscillation  amplitude  should  decrease  when 
the  quenching  is  increased  since  a higher  quenching  rate  would  move  the  lumin- 
osity peak  up  to  v;here  the  electron  density  oscillation  is  less.  This  is  true 


so  long  as  the  peak  luminosity  remains  below  the  F2  peak,  (Porter  and  Tuan,  1973). 

To  understand  why  v/e  should  have  a phase  advance  we  must  consider  the  vari- 
ation in  the  phase  of  the  electron  density  oscillation  with  respect  to  the 
neutral  gravity  wave  as  a function  of  height.  The  solid  curve  in  Fig.  6 shews 
a plot  of  the  theoretical  phase  variation  with  altitude  while  the  dashed  curve 
shows  Thome's  (1968)  observed  results,  (Fig.  6 is  taken  from  Porter  and  Tuan, 
1973).  We  see  that  the  phase  of  the  electron  density  oscillations  advances 


monotonically  with  height.  There  j.s  a region,  however,  below  260  Km  where  an 


increase  in  the  height  of  the  luminosity  will  bring  little  phase  change.  Above 
260  Km,  any  increase  in  the  height  of  the  luminosity  curve  would  bring  a rapid 
advance  in  the  phase  of  the  electron  density  oscillation.  Even  though  Fig.  5 
deals  with  an  integrated  intensity,  since  the  phase  advance  is  mcr.otonic,  an 
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caused  by  gravity  vruves  from  those  caused  by  incident  particle  flux,  if  a 


simultaneous  measurement  cf  total  electron  content  and  airglc.r  intensity’  can 
be  ire  da.  If  no  n re  duct ion  i.  e-chanism  ir  rresont  and  the  disturbance  is  caused 


entirely  by  a gravity  v:a\  c- , then  (1)  there  is  relatively  little  fluctuation 
in  the  total  electron  content  or  in  any  experimental  parameter  that  depends 
only  cn  the  total  number  of  columnar  electrons  rather  than,  say,  how  the  electro 
are  distributed;  (2)  the  rate  of  change  of  total  electron  content  should  remain 
negative.  This  means  that  for  night- time  conditions  and  a gravity  wave  length 
(of  the  order  cf  IbOO  Km  say)  long  compared  with  the  thickness  of  the  F- layer 
so  that  we  may  neglect  horizontal  transport  arising  from  concentration  gradients 


in  the  horizontal  direction,  the  total  columnar  electron  content  should  decree  e 
monoton ically  with  tine.  This  is  easily  understandable  since  the  only  reason 
for  the  little  fluctuations  that  we  do  get  in  the  total  electron  content  is 
due  to  the  variation  in  the  recombination  coefficients  as  the  ionosphere  : .ovc  s 
up  and  down;  (3)  there  is  always  a phase  advance  in  the  airglcw  fluctuations 
with  respect  to  the  rate  of  change  of  total  electron  content.  The  magnitude 
of  the  advance  is  in  general  determined  by  the  quenching.  I’e  should  remark 
here  that  if  a horizontal  time  varying  electrostatic  field  exists  in  the  F- 
region  causing  an  b x B drift  m the  ambient  electrons  only  without  any  non- 
zontal  transport  then  all  three  of  the  above  properties  in  the  behavior  of 
the  air-glow  and  columnar  electron  density  should  also  hold.  They  would  no 
longer  be  true,  however,  if  somehow  additional  electrons  arc  introduced  either 
through  transport  or  some  sort  of  production  mechanism  For  instance , if  the 
disturbance  is  caused  by  an  incident  particle  flux  which  produces  additional 
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ion-electron  pairs,  then  the  rate  of  change  of  total  electron  content  can 
go  positive  since  the  total  columnar  electron  content  can  no  longer  vary 


i:.cr.ct  sntcaxc 


.r.eial.  the  flu ct action  can  aloe  he  far  more 


drastic  (Moschandreas  and  Tuan,  1973). 

As  an  example,  we  consider  the  data  of  Brown  and  Steiger  (1972)  who  have 
simultaneously  measured  the  63C0  a airglcw  and  the  total  electron  content  as 
a function  of  tine  at  Hawaii.  They  found  an  oscillation  in  the  airglow  which 
correlates  with  oscillations  in  the  total  electron  content.  The  negative  tine 
derivative  of  columnar  electron  content  trails  the  air-glow  oscillation  fcy  a 
time  delay  of  about  25  minutes.  We  would  have  at  first  some  reason  to  feel 
that  this  particular  disturbance  could  be  caused  by  a wave  since  (1)  the 
disturbance  appears  to  be  periodic  with  a period  of  about  3 1/2  hours;  (2) 
the  airglow  fluctuation  shows  a phase  advance  of  about  20  to  30  minutes  with 
respect  to  the  rate  of  change  of  total  electron  content.  However,  it  is 
rather  unlikely  that  the  disturbance  is  a wave  since  the  period  is  of  the  order 
of  3 1/2  hours  and  the  disturbance  takes  place  around  the  middle  of  the  night 
when  we  can  forget  about  photo ionisation.  The  disturbance  also  does  not  vary 
monotonicaliy  with  time  and  shows  large  oscillation  amplitudes.  We  may  add 
that  condition  (1)  mentioned  above  can  be  produced  by  any  periodic  disturbance 
and  (2)  can  Le  produced  by  either  a gravity  wave  or  a vertical  drift  without 
any  introduction  of  additional  charged  particles.  There  is,  of  course,  a 
possibility  that  we  may  simultaneously  have  an  electromagnetic  drift  together 
with  an  introduction  of  additional  electrons.  In  fact,  such  an  event  would  be 
quite  consistent  with  both  the  experimental  observations  as  well  as  our  criteria. 
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Conclusion 


V'e  have  examined  the  effect  of  a gravity  wave  in  the  F-region  on  the  electron 
density  and  the  emission  intensity  of  two  forbidden  transitions  whose  mechanism: 


involve  recombination  with  electrons.  Because  of  such  effects  es  quenching  and 


the  lack  of  available  0o+  and  KO  molecules  at  high  altitudes,  the  verticlo  move- 


ment of  the  luminosity  profiles  in  response  to  the  passage  of  the  gravity  wave 
is  much  smaller  than  that  of  the  electron  density  profile.  The  amplitude  of 
oscillation  of  the  airglow,  on  the  other  hand,  is  greater  than  that  of  the  total 
elect  on  density  which  varies  monotcnically  with  time.  Owing  to  the  presence  oi 
a !; agnatic  field  as  well  as  the  shape  of  the  electron  density  profile,  phase 
differences  in  behavior  occur  between  different  heights.  This  effect,  combined 
with  thi  hoi  ;ht  differentiating  effect  cf  quenching  produces  phase  differences 
also  between  airglow  and  ionospheric  parameters.  A combination  of  airglow  and 
ionospheric  measurements  can  be  used  to  differentiate  between  gravity  waves  < 
the  effects  of  incident  particle  flux. 

The  theory  has  been  tested  by  assuming  a gravity  wave  effect  for  some  pub- 
lished data  on  simultaneous  measurements  of  airglow  and  ionospheric  par.-  • ters  . 

The  fit  was  satisfactory. 

This  study  represents  the  first  theoretical  treatment,  as  far  as  we  are  aware, 
of  the  effects  of  gravity  waves  on  airglow.  Y.Te  have  tried  here  to  establish  the 
basic  characteristics  to  be  expected  for  the  luminosity  and  electron  density  pro- 
files under  the  influence  of  a gravity  wave,  and  to  show  how  these  could  be  used 
in  some  canes  to  provide  information  on  processes.  Expansion  of  the  techniques 
should  be  useful  as  a tool  for  the  study  of  both  gravity  waves  and  upper  atmos- 
phere parameters  and  processes. 
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Anrendix 


A good  analytical  approximation  for  the  decay  constant  expressed  as 

a function  of  the  scale  height  factor  p is  given  by: 


.3  D 

, , , v o o 

V*>  = TT 


1 + --i-  ( <*  ^ y(p  + 2g)  - v(l  + h'h 


(A  1) 


A.  derivation  for  this  is  given  by  Tuan  (1968).  Another  totally  different 
derivation  has  been  given  by  Porter  and  Tuan  (1S74).  In  this  expression 
decreases  mcnotonically  v.dth  increasing  p,  where  !<_  P<_  2.  Thus,  the  decay  of 
the  F-lsyer  is  in  general  faster  when  the  gases  are  perfectly  mixed  (P  = 1). 

The  expression  3(q)  in  equation  (2)  is  given  by: 

Aa^  sin©  cos©  y (1  - q) 

B(q)  r; 

2irn  aiH  (A  2) 

As  already  mentioned,  q varies  between  0 and  1/2.  If  q is  set  to  zero  the 
velocity  amplitude  does  not  vary  with  height.  Otherwise,  it  increases  exponent- 
ially. V'e  sco  that  B(q)  vanishes  at  the  magnetic  pole  and  equator  where  we  do 
not  expect  a horizontally  ducted  gravity  wave  to  product  filiation:  in  airj 


The  functions  P..(p,ui,t)  and  P-p(p,u,t)  are  defined  (Porter  and  Tuan,  1974)  by 
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VJh.ile  both  ?_,  and  are  oscillating  functions  of  tine  with  frequency 
\ ~'aY  change  sign  but  ?n  regains  always  positive.  It  is  easy  to  see 


u ->  0,  both  Pn  and  ?,  become  r.onotonic  functions  of  tine  and  all  cscil 
vanish  fron  the  expression  given  by  equation  (2). 

The  functions  F^(p)  and  Cr,  (o)  are  given  by: 
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The  derivation  of  these  function:-  arc  given  in  Tur.n  (1?5S)  and  Pci 
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-hand  side  of  equation  (2) 

to  only  a few  terras  for  a gas  either  perfectly  nixed  or  in  diffusive  o 
and  for  a constant  velocity  amplitude  (a  = C). 
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Figure  Legends 

Tig.  1.  The  comparison  of  the  motion  of  the  electron  density  profile  with 
6300  a luminosity  profile  under  the  influence  of  a gravity  wave. 

Fig.  2.  Time  varieties  c-f  the  0300  ?.  zenith  sir  tier  i®  comparc-d  with  the  tote-1. 

columnar  electron  content  and  the  time  derivative  or  the  total  thee. nan 

content . 

Fig.  3.  Time  variation  of  the  5200  K nitrogen  emission  in  a gravity  wave  for 
the  02  quenching  coefficient  of  Hernandez  and  Turtle  (cashed  curve) 
and  SJanr.c"  et.al.  (solid  curve). 

fig.  *4.  The  experimental  observations  of  Luchs  shewing  simultaneous  variations 
in  foF2,  h'F2  and  the  6300  /•.  airglow  are  compared  with  theoretical 
calculations  (solid  curves)  based  or.  the  assumption  of  a gravity  wave 
with  a period  of  approximately  3 hours. 


The  oscillations  of  the  6500  i.  line  for  two  nitrogen  quenching 
coefficients,  showing  how  the  larger  quenching  coefficic.  x (dashed 
curve)  reduces  oscillation  ana  introduces  a phase  advance. 


Comparison  of  the  theoretical  calculations  of  the  variation  with 
altitude  in  the  phase  of  the  ionospheric  oscillations  with  Thome’s 
(1268)  observed  values. 
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Although  the  experimental  observations  of  airglow  intensity  fluctuations  and  their  possible  connections 
with  traveling  ionospheric  disturbances  have  been  made  by  many  investigators,  very  little  attention  has 
been  given  to  a theoretical  analysis  of  such  ionospheric  disturbances  on  airglow.  In  this  paper  we 
specifically  examine  the  effect  of  an  atmospheric  gravity  wave  on  the  6300-A  O 1 and  the  5200-A  N I emis- 
sion in  mid-latitudes.  We  will  also  show  that  for  the  6300-A  O I if  we  assume  reasonable  values  for  the 
ionospheric  parameters  as  well  as  a gravity  wave  with  a period  of  approximately  3 h.  the  resultant  fluc- 
tuations in  airglow  are  consistent  with  Dachs's  observed  data  as  well  as  with  the  simultaneously  measured 
variations  in  the  ionospheric  parameters. 


An  association  of  airglow  intensity  variations  with  traveling 
ionospheric  disturbances  has  been  established  by  Barbier 
[1964,  1965]  for  the  6300-A  O I emission  and  by  Weill  and 
Christophe-Glaume  [1967]  for  both  the  6300-A  O I and  the 
5200-A  N I emission  by  a diicct  comparison  of  the  two  types 
of  measurements.  In  addition,  some  published  data  [Silver- 
man.  1962;  Okuda.  1962;  Dachs.  1968]  can  in  retrospect  be  seen 
as  having  been  probably  due  to  similar  causes.  The  data  of 
Dachs  include  observations  on  both  5577-A  O I and  6300-A  O 
I as  well  as  observations  of  the  ionospheric  parameters  f0F2 
and  It'hf  and  are  thus  particularly  valuable  for  analysis  from 
the  point  of  view  of  a possible  wave  origin  of  the  variations. 
Disturbances  in  airglow  can  of  course  arise  from  a variety  of 
causes,  and  simultaneous  measurements  of  other  parameters 
are  needed  to  establish  the  cause  of  a disturbance. 

Over  the  past  few  years  considerable  theoretical  and  ex- 
perimental work  has  been  carried  out  on  the  effects  of  at- 
mospheric gravity  waves  on  ionospheric  parameters  The 
effects  of  these  waves  on  optical  emissions,  however,  have 
received  no  attention  from  a theoretical  viewpoint.  A suitable 
analytic  model  for  such  a treatment  has  been  presented 
recently  by  Porter  and  Tuan  [1974].  This  model  for  the 
response  of  the  F layer  to  gravity  waves  differs  fundamentally 
from  the  usual  theoretical  treatment  in  several  respects.  Of 
these  the  most  pertinent  for  airglow  applications  are  the 
following: 

I.  The  analytic  model  is  fully  time  dependent.  One  result 
of  such  a model  is  that  given  a sinusoidal  neutral  gravity  wave, 
the  response  of  the  electron  density  oscillation  is  no  longer  in 
general  sinusoidal  and  the  time  behavior  varies  with  height 
This  makes  it  more  suitable  for  studying  both  the  general 
decay  of  the  F layer  and  the  airglow  at  night  as  well  as  such 


detailed  behavior  as  differences  in  phase  between  the  os- 
cillations of  electron  density  and  the  airglow. 

2.  The  analy  tic  model  remains  valid  even  if  the  percentage 
oscillations  in  the  electron  density  or  airglow  are  large  in  com- 
parison with  ambient  values.  This  condition  is  especially  im- 
portant lor  the  6300-A  O I emission,  whose  luminosity  profile 
peaks  at  a position  (about  one  scale  height  below  the  F,  peak) 
where  the  electron  density  and  airglow  oscillations  are  rather 
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In  this  paper  we  first  outline  the  computation  for  the  6300-A 
O I emission  intensity  during  the  passage  of  a gravity  wave. 
For  long-period  gravity  waves  the  continuity  equation  from 
which  the  emission  intensity  of  the  6300-A  O 1 can  be 
calculated  reduces  to  the  quasi-equilibrium  theories  of  Peter- 
son el  al.  [1966],  in  which  the  production  of  the  excited  0(7)) 
states  is  assumed  to  equal  the  loss  through  optical  emission 
and  quenching.  Comparison  of  the  theoretically  calculated 
changes  in  both  electron  density  and  luminosity  profiles  will 
serve  to  bring  out  the  essential  features  of  the  theory  . 

We  shall  next  consider  the  time  behavior  of  the  total  elec- 
tron content  as  compared  with  the  airglow  oscillations.  It  will 
be  shown  that  the  total  electron  content  shows  less  oscillatory 
behavior  and  that  only  by  plotting  the  time  derivative  of  the 
total  electron  content  can  the  oscillations  be  more  significant. 
However,  such  fluctuations  arc  detectable  [Davis  and  daRosa. 
1969],  The  reduced  oscillation  in  the  total  electron  content, 
coupled  with  the  large  optical  emission  oscillations,  provides 
us  with  a possible  criterion  for  distinguishing  between  long- 
period  gravity  wave  effects  and  effects  produced  by  other 
sources  such  as  particle  participation. 

We  will  then  consider  the  effect  of  a gravity  wave  on  the 
5200-A  N I emission  intensity  . Here  we  can  no  longer  use  the 
equilibrium  theory  that  we  have  been  using  for  the  6300-A 
line,  since  at  least  at  highei  altitudes  the  pioduction  of  N(7>) 
states  can  no  longer  be  balanced  by  either  the  radiative  loss  or 
the  quenching  We  will  also  consider  the  importance  of  elec- 
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tron  quenching,  which  is  nol  significant  for  the  6300-A  O I 
line,  but.  us  we  shall  see,  is  much  more  important  for  the 
nitrogen  line. 

The  theory  derived  here  will  then  be  used  lor  the  interprela- 
- lion  of  the  experimental  results  of  Dachs  [1968]  at  Tsutneb. 
South-West  Africa,  These  include  simultaneous  measurements 
of  /..F,.  h' and  the  airglow  intensity  variations  with  time, 
which  allow  for  a test  of  both  aspects  of  the  theory. 

Finally,  we  shall  also  study  the  effect  of  varying  the 
quenching  rales  on  the  phase  of  the  oscillations.  We  will  con- 
sider this  phase  difference  in  relation  to  the  effect  of  other 
mechanisms  that  can  also  cause  a phase  difference. 

Time  Behavior  ok  Electron  Density  and 
6200-A  O 1 Luminosity  Profiles 

Theoretical  treatments  of  the  6300-A  O 1 emission  taking  ac- 
count of  different  factors  and  calculating  various  parameters 
are  available  from  a number  of  authors  [e.g.,  Chamberlain, 
1958;  Lagos  el  al..  1963;  Lagos,  1964;  Peterson  el  ai,  1966; 
Tuan,  1969],  The  following  treatment  leans  most  heavily  on 
the  work  of  Peterson  et  al. 

We  make  the  following  assumptions  (for  their  justification 
see  Peterson  et  al.  [1966]):  (I)  the  time  derivative  term  in  the 
continuity  equation  for  the  density  of  oxygen  atoms  in  their 
O I'D)  states  may  be  neglected.  (2)  the  diffusion  term  is 
neglected,  and  (3)  we  assume  chemical  equilibrium  for  the 
production  and  loss  of  02*  molecules  and  neglect  the  02*  and 
NO*  molecules  in  the  equation  for  charge  neutrality.  We  then 
obtain  the  following  well-known  expression  for  the  volume 
emission  rate  <630(l  for  the  6300-A  line: 

/ . _ ( 7 [O,,  ] A', 

8300  ’ ’ V\  A, J 1 + (IM„)(S,[/V2]  + Si  N,) 

where  i ; is  the  efficiency  for  recombination.  7 is  the  rate 
coefficient  for  charge  transfer  between  O and  02' , and  ,V,  and 
S 2 are  the  nitrogen  and  electron  quenching  coefficients,  respec- 
tively. The  02‘  profiles  obtained  by  Schunk  and  Walker  [1973] 
could  later  be  used  in  (I).  However,  since  the  same  chemical 
equilibrium  and  charge  neutrality  approximations  have  been 
used  in  the  expression  for  N,.  [Porter  and  Tuan.  1974],  we  have 
decided  to  maintain  consistency  and  use  (I)  for  the  emission 
rate.  Since  our  primary  interest  is  in  the  basic  dynamic 
response  of  the  airglow  to  gravity  waves  rather  than  any 
detailed  behavior,  we  have  also  neglected  the  contribution 
from  dissociative  recombination  of  NO'  [ Silverman , 1970]  as 
well  as  the  cascading  from  the  '.V  states  in  the  production  of  'P 
states  [ Rishheth  and  Oarriott,  1969],  We  will  also  neglect  the 
time  variation  in  O,  and  N2  caused  by  the  gravity  waves.  In 
reality,  the  equations  for  such  neutral  atmospheric  con- 
stituents are  really  coupled  to  the  equation  for  the  ion-electron 
plasma.  To  simplify  all  these  processes,  we  have  used 
stationary  models  for  02  and  N2.  In  practice  the  ap- 
proximations should  be  good,  since,  as  was  already  men- 
tioned, the  percentage  variation  in  the  neutral  gas  caused  by 
the  gravity  wave  is  usually  very  much  less  than  the  percentage 
variation  for  the  ion-electron  plasma,  especially  in  important 
regions  such  as  the  peak  of  the  luminosity  profile  By  far  the 
most  important  contribution  to  the  oscillation  in  the  airglow 
arises  from  changes  in  the  electron  density  \,  in  (I) 

Most  recent  theoretical  treatments  for  the  time-dependent 
behavior  of  / layer  electron  density  (e  g.,  Testud  and  Francois. 
1971:  klnstcrmevcr.  |972«.  h\  under  the  influence  of  gravity 
waves  arc  essentially  based  on  perturbation  theory  in  which 


the  unperturbed  ionosphere  is  assumed  to  be  a stationary 
L hapntai.  function  independent  of  time.  Also,  in  perturbation 
treatments  the  time  dependence  for  the  plasma  is  assumed  to 
be  sinusoidal  if  the  gravity  wave  is  sinusoidal,  and  a simple 
product  form  for  the  time  and  spatial  dependence  is  usually 
assumed. 

The  luminosity  peak  for  the  nighttime  6300-A  line  occurs 
about  one  scale  height  below  the  h2  peak  [Tuan.  1969],  where 
the  electron  density  oscillation  can  be  larger  [Porter  and  Tuan. 
1974],  As  an  example,  for  the  gravity  wave  that  Thome  [ 1 968 ] 
analyzed  the  oscillation  in  electron  density  at  this  height  would 
be  between  20  and  25  H of  the  ambient  value.  If  we  go  down  by 
another,  say,  three  quarters  of  a scale  height,  where  the 
luminosity  is  still  very  considerable,  the  oscillation  in  electron 
density  rises  to  about  45-50<!»  of  its  ambient  value;  thus  any 
perturbation  approach  is  invalidated. 

An  analytic  model  for  the  electron  density  A ,(z,  t ) under  the 
influence  of  a gravity  wave  has  been  developed  [Porter  and 
Titan.  1974]  that  seems  appropriate  for  analyzing  airglow  and 
is  valid  for  large  disturbances.  The  results  are  fully  time  depen- 
dent and  may  be  used  to  study  the  general  decay  (including 
changes  in  shape)  of  both  the  F layer  and  the  airglow  in  the 
absence  of  production  of  electrons.  Finally,  we  can  use  this 
model  to  study  the  time  variation  of  the  height  of  the  h\  and 
luminosity  peaks,  a study  that  would  be  difficult  in  a partially 
time-dependent  theory  that  does  not  allow  for  any  changes  in 
shape  in  the  ambient  ionosphere.  In  this  model  we  have 
neglected  time  variations  in  the  temperature,  diffusion,  and 
recombination  coefficients  (all  of  which  can  be  incorporated 
without  too  much  difficulty),  since  it  can  be  shown  that  at  least 
for  long-period  gravity  waves  [Klostermever.  19726:  Porter  and 
Tuan.  1974]  the  percentage  variation  of  these  parameters  with 
time  is  small  over  most  parts  of  the  F layer.  For  the  case  of 
long-period  horizontally  ducted  gravity  waves,  the  expression 
for  the  elctron  density  A ,.(r.  1),  [Porter  and  Tuan.  1974.  equa- 
tion (7)],  which  fully  includes  effects  due  to  diffusion  and 
recombination,  is  given  by 

"'=C  exP  H/,y+  2 cxP  (-//)]} 

exp  [—  X„(r  — r„)  + B(q)  sin  u(t  — 7,)] 

— 5Z  p-(p,  0>.  r)Fn(p) 

1 L 

X L.  "[»  «n  ( -„)]}  (!) 

where  It  h„  is  the  difference  between  the  actual  height  It 
and  some  arbitrary  reference  height  6„.  H is  the  scale  height;  /> 
is  the  scale  height  factor;  \„(/>)  is  a decay  constant  (sec 
appendix);  u>  is  the  angular  frequency  of  the  gravity  wave,  q is  a 
parameter  (0  q ■ '.•)  that  controls  the  damping  of  the 

gravity  wave  (when  q we  obtain  the  undamped  Hun  s 

[I960]  mode'll.  IHq l is  ,1  coefficient  depending  on  gravitv  wave 
and  ionospheric  parameters  (see  appendix  i;  I,  is  some  initial 
time  giving  us  the  phase  ol  the  gravity  wave  al  1 i„;  n 

21  ltd,,  /)„)'•'.  where  A,  and  />  arc  the  recombination  and 
diffusion  coefficients  at  h„.  f)  is  the  angle  between  the  magnetic 
field  and  the  zenith,  I is  the  amplitude  of  the  gravity  wave  at 
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/'„(/>.  w.  i)  and  Pip,  «.  I)  arc  oscillating  functions  of  lime 
defined  in  the  appendix;  /■„(/>)  and  (<„(</)  are  factors  that  con- 
trol the  dependence  on  the  scale  height  factor  and  the  damping 
of  the  gravity  wave,  respectively,  and  are  given  in  the  appen- 
dix; L„  ' 2 is  a Laguerre  polynomial  of  half-integral  order. 

We  see  at  once  from  (2)  that  if  we  ignore  the  contribution 
from  the  infinite  sum  on  the  right-hand  side,  we  have  a Chap- 
man function  multiplied  by  an  exponential  damping  together 
with  oscillations.  This  is  as  would  be  expected,  since  for  the 
nighttime  situation  considered  here,  production  can  be  ex- 
cluded. (We  may  easily  include  production  if  it  is  required.) 

The  infinite  sum  on  the  right-hand  side  of  (2)  reduces  to  a 
few  terms  for  a gas  that  is  either  perfectly  mixed  or  is  in 
diffusive  equilibrium  and  for  a constant  velocity  amplitude 
(independent  of  height)  for  the  gravity  wave  (</  = 0)  (see 
appendix). 

In  the  case  when  there  is  no  gravity  wave  (A  = 0)  or  when  we 
are  at  the  geomagnetic  pole  (II  = 0)  or  equator  (6  = jt/2).  B(q) 
- 0 (see  appendix),  and  ( I ) reduces  to  the  expression  for  an  un- 
disturbed ionosphere  [Moschandreas  and  Tuan.  1974].  The 
decay  constant  \„(p ) describes  the  general  Flayer  decay  in  the 
absence  of  disturbance  and  w as  originally  derived  through  per- 
turbation theory  | Tuan.  1968],  In  general,  \„(p)  increases  when 
the  gases  are  in  diffusive  equilibrium  (/>  = 2). 

The  expression  given  by  (2)  provides  the  correct  phase 
behavior  as  well  as  the  F2  peak  motion,  as  is  seen  from  a com- 
parison with  Thome's  [1968]  observations  [Purler  and  Tuan. 
1974],  By  inserting  the  above  expression  for  ,V,.(r.  /)  in  (I)  we 
obtain  the  luminosity  profiles  for  the  6300-A  O 1 shown  in 
figure  I.  Here  we  have  adopted  ionospheric  and  gravity  wave 
parameters  similar  to  those  used  to  til  Thome's  [1968]  obser- 
vations except  for  the  scale  height  factot  for  which  we  have 
now  used  p I 75  in  order  to  bring  the  gases  closer  to  diffusive 
equilibrium  as  one  might  expect  for  the  I-  region.  We  shall 
assume  that  A I I 10  s 1 [ Chamberlain , 1961],  S,  5 X 
10  " cm3  s 1 | Peterson  and  Van  /.and!,  1969],  and  S2  1.7  X 
10  " cm1  s 1 [Seaton . 1 956J . 

In  figure  I a comparison  of  the  electron  density  profile  dur- 
ing the  passage  of  a gravity  wave  with  the  movements  of  the 
luminosity  profile  seems  to  show  that  (I)  in  general,  the 
luminosity  profile  peaks  about  one  scale  height  below  the  F2 
peak  as  it  does  in  the  undisturbed  time-dependent  theory 
[ Tuan.  1969],  (2)  I hour  later,  as  the  I , peak  moves  down  by  20 
km  losing  all  the  time  in  electron  density  magnitude,  the 
luminosity  peak  remains  relatively  unchanged  in  position  but 
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greatly  enhanced  in  intensity  (figure  I.  center),  and  (3) 
another  hour  later,  as  the  Ft  peak  moves  up  again,  the  inten- 
sity of  peak  luminosity  again  drops,  there  being  little  rise  in 
position  (figure  I,  right).  If  we  had  plotted  the  luminosity 
profiles  withot  uenching,  the  peak  luminosity  would  have 
moved  a consioerable  distance  (about  15  km)  downward  when 
the  T2  peak  came  down.  This  conclusion  implies  that  the 
luminosity  peak  is  essentially  kept  up  by  the  rapid  increase  in 
quenching  at  lower  altitudes.  The  reason  that  it  would  only 
move  15  km  rat:  r than  the  20  km  for  the  F2  peak  is  the  rapid 
drop  in  available  t )./  for  the  production  of  'D  states  at  higher 
altitudes,  which  lends  to  keep  the  luminosity  peak  down.  Thus 
the  combined  effect  of  both  is  to  keep  the  luminosity  peak 
from  having  larger  oscillations  in  the  vertical  direction. 

Oscillations  in  Totai  Electron 
Content  and  6300-A  O I 

I o obtain  th  . tal  zenith  airglow  intensity  for  the  6300-A  O 
I.  we  integrate  (I)  over  height.  The  results  are  plotted  against 
local  time  ie  I igtu  e 2.  We  have  also  plotted  the  total  columnar 
electron  content  a . well  as  the  negative  time  deriv  ative  of  this 
total  electron  content.  The  parameters  used  here  are  the  same 
as  those  adopted  for  f igure  1.  One  general  feature  is  that  all 
three  curves  decay  with  time,  an  expected  result,  since  no  ion 
pair  production  mechanisms  have  been  included.  The  most 
significant  feature  is  that  although  there  is  considerable  oscilla- 
tion with  time  in  the  airglow  intensity,  the  total  electron  con- 
tent shows  ' Intis  ly  little  fluctuation.  The  lack  of  fluctuation 
in  the  lot  i1  n content  is  readily  understandable.  We  are 

lii si  •..• . v, ah  a nighttime  situation,  when  wc  may 

neglect  ph  iclion.  The  fluctuation  , arise  when  the  tons 

are  moved  il  itely  into  rarified  atmosphere,  where  the 

recombinat-  less,  and  dense:  atmosphere,  where  the 

recombination  treater.  The  fluctuation  i further  inhibited 
by  phase  cancellation  between  the  motion  of  ions  above  and 
below  the  / peak  Since  we  are  dealing  with  a horizontally 
dueled  neutral  wave  with  no  vertical  ph  ixc  progression,  the 
phase  variation  of  ions  depends  almost  entirely  on  the  shape  of 
the  F layer,  f igure  6 shows  how  the  phase  of  the  ion  motion 
above  the  vicinity  of  the  F2  peak  is  in  general  opposite  to  that 
below  the  peak  We  have  here  an  example  of  possible  ph  e 
cancellation  even  when  the  gravity  wave  does  not  possess 
phase  variation  along  the  line  of  sight  <■  ot 
] 1970]  have  show  i how  phase  cancellation  mav  v o Ir 
phase  vari  i the  gravity  wave  ilscM 
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irons  over  . and  Topping  the  pr  hi 
result  that  the  columnar  recombm  it 
negative  of  the  tale  of  ch.ingi 
since  the  divergence  term  in  th 
from  Ciauss  theorem  I hus  i - f 
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6300  A ZENITH  AIRGLOW 


Fig.  2.  Time  variation  of  the  6300-A  zenith  airglow  compared  with 
the  total  columnar  electron  content  and  the  time  derivative  of  the  total 
electron  content. 

quenching  as  we  have  just  mentioned  or  to  introduce  some 
production  process.  This  point  will  be  discussed  further  later. 

5200-A  Nitrogen  Line 

The  5200-A  N I is  produced  in  the  F region  by  dissociative 
recombination  in  a manner  similar  to  the  production  of  the 
6300-A  line.  Analysis  of  the  effects  of  a gravity  wave  on  this 
line  is  of  special  interest  because  of  the  very  long  lifetime  of  the 
nitrogen  1DS  2 states.  We  can  no  longer  use  the  quasi- 
equilibrium theories,  especially  at  higher  altitudes,  where 
quenching  can  no  longer  compensate  for  the  very  slow 
radiative  loss.  If  we  neglect  diffusion  [ Hernandez  and  Turtle, 

1969],  the  continuity  equation  for  the  atomic  nitrogen  density 
in  the  ’£>  state  N(’D),  is  given  by 

= P - ( A, „+  K,[02]  + KtN.)NCD)  (3) 

In  (3),  P is  the  production  rate  for  N(’£>),  K,  and  Kt  are  the 
quenching  coefficients  for  02  and  electrons,  respectively,  A ^ is 
the  Einstein  coefficient  for  N(’D).  If  we  go  down  to  low 
enough  altitudes,  where  the  quenching  is  sufficiently  strong  to 
balance  the  production  rate,  we  may  neglect  the  time 
derivative  term,  and  the  theory  reduces  to  a quasi-equilibrium 
theory  exactly  analogous  to  that  of  Peterson  et  at.  [1966).  For 
the  production  rate  we  again  assume  that  the  production  of 
N(’/3)  through  dissociative  recombination  of  NO*  is  balanced 
by  charge  transfer  between  O’  and  N2.  In  this  way,  the 
production  rate  of  N(’£>)  states  varies  with  height  as  the 
product  of  neutral  N2  concentration  and  N,. 

The  loss  rate  of  N(2/))  is  represented  by  the  three  terms  on 
the  right-hand  side  of  (3),  corresponding  to  radiative  loss,  02 
quenching,  and  quenching  due  to  superelastic  collisions  with 
the  electrons.  We  have  assumed  that  A ^ = 1/26  h 1 
[ Chamberlain , 1961]  and  K,  = 8 X I0'1”  cm’  s“‘  \Seaton, 

1956].  For  the  rate  coefficient  K,  with  respect  to  02  quenching 
we  have  tried  both  K,  = 1.4  X 10  " cm’ s 1 [Slanger  et  al., 

1971]  and  #£,  = 1.4  X 10  ” cm*  s ' ( Hernandez  and  Turtle, 

1969],  The  two  values  differ  by  an  order  of  magnitude,  and  * 
there  is  reason  to  believe  that  the  value  given  by  Hernandez 
and  Turtle  may  be  on  the  small  side  (G.  Hernandez,  private 
communication,  1973).  Equation  (3)  may  be  simply  in- 
tegrated. subject  to  the  initial  boundary  conditions.  From  the 
N (’/))  concentration  we  may  readily  calculate  the  columnar 
zenith  airglow  intensity  for  the  (*S,  2-’/7,  forbidden  op- 
tical transition  as  a function  of  time.  The  results  are  shown  by 


the  solid  curve  in  Figure  3 for  the  quenching  coefficient  of 
Slanger  et  al  [1971]  and  by  the  dotted  curve  for  the  coefficient 
used  hy  Hernandez  and  Turtle  [1969].  For  both  curves  we  have 
used  a scale  height  factor  of  p = 1.75.  All  other  gravity  wave 
parameters  are  the  same  as  those  used  to  fit  Thome's  results 
[Porter  and  Tuan,  1974], 

The  three  important  features  are  that  (I)  the  oscillations  for 
the  higher  02  quenching  rate  of  Slanger  el  at.  [1971]  are  less 
pronounced  than  those  produced  by  the  lower  quenching 
coefficient  of  Hernandez  and  Turtle  [1969],  (2)  there  is  a phase 
advance  with  the  higher  quenching  coefficient  of  between  10 
and  20  min  (the  reasons  for  this  will  be  given  in  a later  section), 
and  (3)  the  luminosity  peak  for  the  higher  quenching 
coefficient  is  in  general  about  20  km  higher  than  the  peak  for 
the  lower  coefficient. 

In  comparing  these  results  with  the  behavior  of  the  6300-A 
line,  there  are  the  following  important  differences:  The  elec- 
tron quenching  for  the  nitrogen  line  is,  relatively  speaking, 
much  more  important  than  that  for  the  oxygen  line.  If  we  use 
the  coefficient  of  Slanger  et  al  [1971],  the  quenching  due  to 
molecules  for  the  nitrogen  line  is  down  by  more  than  an  order 
of  magnitude  from  the  quenching  for  oxygen,  owing  to  the 
combined  effect  of  a smaller  quenching  coefficient  and  a 
smaller  abundance  of  02  molecules  than  of  N2  molecules.  The 
difference  in  electron  quenching  is  only  about  a factor  of  2. 

The  relatively  more  important  electron  quenching  has  the 
effect  of  holding  down  the  5200-A  luminosity  peak,  since  the 
available  electrons  for  quenching  increase  until  the  F 2 peak, 
well  above  the  luminosity  peaks  of  both  the  6300-A  and  the 
5200-A  line. 

The  net  result  is  that  although  one  may  at  first  expect  the 
5200-A  luminosity  peak  to  be  well  above  the  oxygen  peak,  ow- 
ing to  a much  longer  half  life,  the  reduced  molecular 
quenching  at  lower  altitudes  coupled  with  a relative  increase  in 
electron  quenching  (unimportant  for  the  6300-A  line)  at  higher 
altitudes  would  bring  down  the  peak  position  for  the  5200-A 
line  to  height  levels  comparable  to  those  for  the  6300-A  line. 
This  is  not  inconsistent  with  the  observations  of  Weill  and 
Christophe-Giaume  [1967],  who  argued  from  geometric  con- 
siderations that  the  5200-A  peak  is  about  12  km  above  the 
6300-A  peak  if  the  speed  of  the  traveling  ionospheric  distur- 
bance is  assumed  to  be  the  same  at  the  peaks  of  both 
emissions.  We  may  conclude  that  since  both  peaks  occur  at 
nearly  the  same  height,  the  oscillation  amplitudes  should  be 
comparable.  This  conclusion  is  borne  out  by  comparing 
Figure  3 with  Figure  2. 


Fig  3.  Time  variation  of  the  5200-A  nitrogen  emission  in  a gravity 
wave  for  the  O,  quenching  coefficient  of  Hernandez  and  Turtle  |I969| 
(dashed  curve)  and  Slanger  el  al  |I97I|  (solid  curve) 
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Comparison  With  Simultaneous  Observations  of 
JoF2,  h'F2,  and  6300-A  O I Emission 

We  have  mentioned  earlier  that  some  published  data  seem 
likely  to  have  exhibited  the  effects  of  gravity  waves.  In  this  sec- 
tion we  use  some  of  Dachs' s [1968]  simultaneous  measurement 
of  ionospheric  parameters  and  the  6300-A  O I emission  to  il- 
lustrate the  theory  by  fitting  the  observations  to  an  assumed 
gravity  wave. 

Dachs  [1968]  has  made  intensity  measurements  of  the  6300- 
and  5577-A  emission  lines  as  a function  of  time.  These 
measurements  together  with  simultaneous  measurements  of 
l„F2  and  h'F2  seem  to  show  considerable  correlation,  especially 
between  the  6300-A  line  and  the  ionospheric  parameters.  In 
Figure  4 we  have  shown  a comparison  between  our  theoretical 
results  and  the  experimental  values  for  f0Ft,  h'F2 , and  the  6300- 
A airglow.  For  the  theoretical  parameters  we  have  chosen  p = 
1.5,  H = 65  km,  the  recombination  coefficients  0(200  km)  = 
0.34  h ',  and  the  diffusion  coefficient  D(200  km)  = 6.1  X 103 
km2  h ‘.  These  are  somewhat  different  from  the  parameters 
chosen  to  fit  Thome's  observation,  but  considering  the  fact 
that  we  have  completely  ignored  production  processes  (the 
observations  are  made  between  1800  and  2400  LT),  they  are 
the  right  order  of  magnitude.  For  the  gravity  wave  we  assume 
the  velocity  amplitude  A( 200  km)  = 35.2  m s“‘  and  the  damp- 
ing factor  q = 0.49,  which  is  rather  close  to  Hines’s  [I960]  un- 


THEORY 


Fig  4 Experimental  observations  (open  circles)  of  Dachs  showing 
simultaneous  variations  in  f„F„  h'F,,  and  the  6300-A  airglow  com- 
pared with  theoretical  calculations  (solid  curves)  based  on  the  assump- 
tion of  a gravity  wave  with  a period  of  approximately  3 h. 


Fig.  5.  Oscillations  of  the  6300-A  line  for  two  nitrogen  quenching 
coefficients,  showing  how  the  larger  quenching  coefficient  (dashed 
curve)  reduces  oscillation  and  introduces  a phase  advance. 

damped  model.  The  quenching  at  200  km  is  assumed  to  be 
0.24  s ',  and  the  electron  quenching  is  neglected.  The  quantum 
efficiency  for  emission  from  the  'D  is  assumed  to  be  0.5. 

In  view  of  the  rather  good  agreement  between  theory  and 
experiment  for  both  the  airglow  and  the  ionospheric 
parameters  one  can  readily  conclude  that  even  without  a total 
electron  content  measurement,  we  have  a gravity  wave  passing 
through  with  a period  of  between  2 and  3 h.  The  comparison 
between  the  airglow  intensity  and  the  h'F2  shows  all  the  ex- 
pected phase  correlations;  i.e.,  when  the  airglow  intensity  is  at 
a minimum,  the  F2  peak  is  approximately  at  a maximum. 

Effect  of  Quenching  and  the  Criteria  for 
Different  Disturbances 

As  was  already  mentioned  in  a previous  section,  quenching 
can  cause  a phase  difference  in  the  oscillation  of  the  airglow  in- 
tensity. Figure  5 shows  the  6300-A  O I for  two  N2  quenching 
coefficients.  The  solid  curve  uses  a quenching  coefficient  S,  = 
1 .0  X 10  11  cm3  s ',  whereas  the  dashed  curve  uses  S,  = 2.7  X 
10  " cm3  s All  other  parameters  are  the  same  as  those  used 
to  fit  Dachs' s [1968]  experimental  observations.  The  highei 
quenching  coefficient  seems  to  produce  fewer  oscillations  as 
well  as  a phase  advance  by  about  half  an  hour.  It  is  easy  to  un- 
derstand why  the  oscillation  should  decrease  when  the 
quenching  is  increased,  since  a higher  quenching  rate  would 
move  the  luminosity  peak  up  to  a point  where  the  electron 
density  oscillation  is  less.  This  is  true  so  long  as  the  peak 
luminosity  remains  below  the  F2  peak  [Porter  and  Tuan.  1974], 

To  understand  why  we  should  have  a phase  advance,  we  must 
consider  the  variation  in  the  phase  of  the  electron  density  os- 
cillation with  respect  to  the  neutral  gravity  wave  as  a function 
of  height.  The  solid  curve  in  Figure  6 shows  a plot  of  the 
theoretical  phase  variation  with  altitude,  and  the  dashed  curve 
shows  Thome's  [1968]  observed  results.  (Figure  6 is  taken  from 
Porter  and  Tuan  [1974].)  We  see  that  the  phase  of  the  electron 
density  oscillations  advances  monotonically  with  height.  There 
is  a region,  however,  below  260  km  where  an  increase  in  the 
height  of  the  luminosity  will  bring  little  phase  change.  Above 
260  km  any  increase  in  the  height  of  the  luminosity  curve 
would  bring  a rapid  advance  in  the  phase  of  the  electron  den- 
sity oscillation.  Even  though  Figure  5 deals  with  an  integrated 
intensity,  since  the  phase  advance  is  monotonic,  an  upward 
shift  of  the  entire  luminosity  profile  must  always  produce  a 
phase  advance. 

The  analysis  seems  to  suggest  a means  for  distinguishing  air- 
glow fluctuations  caused  by  gravity  waves  from  those  caused 
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Fig.  6.  Theoretical  calculations  of  the  variation  with  altitude  in  the 
phase  of  the  ionospheric  oscillations  compared  with  Thome’s  (1968) 
observed  values. 


by  incident  particle  flux,  if  a simultaneous  measurement  of 
total  electron  content  and  airglow  intensity  can  be  made.  If  no 
production  mechanism  is  present  and  the  disturbance  is  caused 
entirely  by  a gravity  wave,  then  the  following  observations  are 
true: 

1.  There  is  relatively  little  fluctuation  in  the  total  electron 
content  or  in  any  experimental  parameter  that  depends  only 
on  the  total  number  of  columnar  electrons  rather  than,  say,  on 
how  the  electrons  are  distributed. 

2.  The  rate  of  change  of  total  electron  content  should  re- 
main negative.  Thus  for  nighttime  conditions  and  a gravity 
wavelength  that  is  sufficiently  long  (of  the  order  of  1500  km, 
say)  in  comparison  with  the  thickness  of  the  F layer  that  we 
may  neglect  horizontal  transport  arising  from  concentration 
gradients  in  the  horizontal  direction,  the  total  columnar  elec- 
tron content  should  decrease  monotonically  with  time.  This 
result  is  easily  understandable,  since  the  only  reason  for  the  lit- 
tle fluctuations  that  we  do  get  in  the  total  electron  content  is 
the  variation  in  the  recombination  coefficients  as  the 
ionosphere  moves  up  and  down. 

3.  There  is  always  a phase  advance  in  the  airglow  fluc- 
tuations with  respect  to  the  rate  of  change  of  total  electron 
content.  The  magnitude  of  the  advance  is  in  general  deter- 
mined by  the  quenching. 

We  should  remark  here  that  if  a horizontal  time-varying 
electrostatic  field  exists  in  the  F region  causing  an  E x B drift 
in  the  ambient  electrons  only,  without  any  horizontal 
transport,  then  all  three  of  the  above  properties  in  the  behavior 
of  the  airglo u and  columnar  electron  density  should  also  hold. 
They  would  no  longer  be  true,  however,  if  somehow  additional 
electrons  were  introduced  either  through  transport  or  some 
sort  of  production  mechanism.  For  instance,  if  the  disturbance 
is  caused  by  an  incident  particle  flux  that  produces  additional 
ion-electron  pairs,  then  the  rate  of  change  of  total  electron 
content  can  go  positive,  since  the  total  columnar  electron  con- 
tent can  no  longer  vary  monotonically  with  time.  In  general, 
the  fluctuation  can  also  be  far  more  drastic  [ Moschandreas  and 
Tuan.  1974). 

As  an  example  we  consider  the  data  of  Brown  and  Steiger 
(I972).  who  have  simultaneously  measured  the  6300-A  airglow 
and  the  total  electron  content  as  a function  of  time  at  Hawaii. 
They  found  an  oscillation  in  the  airglow  that  correlates  with 
oscillations  in  the  total  electron  content.  The  negative  time 
derivative  of  columnar  electron  content  trails  the  airglow  os- 


cillation by  a time  delay  of  about  25  min.  We  have  at  first  some 
reason  to  fee!  that  this  particular  disturbance  could  be  caused 
by  a wave,  since  ( 1 ) the  disturbance  appears  to  be  periodic  with 
a period  of  about  31 : h,  and  (2)  the  airglow  fluctuation  shows  a 
phase  advance  of  about  20-30  min  with  respect  to  the  rale  of 
change  of  total  electron  content.  However,  it  is  rather  unlikely 
that  the  disturbance  is  a wave,  since  the  period  is  of  the  order 
of  31:  h and  the  disturbance  takes  place  around  the  middle  of 
the  night,  when  we  can  forget  about  photoionization.  The  dis- 
turbance also  does  not  vary  monotonically  with  time,  and  it 
shows  large  oscillation  amplitudes.  We  may  add  that  the  first 
condition  mentioned  above  can  be  produced  by  any  periodic 
disturbance,  and  the  second  condition  can  be  produced  by 
either  a gravity  wave  or  a vertical  drift  without  any  introduc- 
tion of  additional  charged  particles.  There  is  of  course  a 
possibility  that  we  may  simultaneously  have  an  elec- 
tromagnetic drift  together  with  an  introduction  of  additional 
electrons.  In  fact,  such  an  event  would  be  quite  consistent  with 
both  the  experimental  observations  and  our  criteria. 

Conclusion 

We  have  examined  the  effect  of  a gravity  wave  in  the  F 
region  on  the  electron  density  and  on  the  emission  intensity  of 
two  forbidden  transitions  whose  mechanisms  involve  recom- 
bination with  electrons.  Because  of  such  effects  as  quenching 
and  the  lack  of  available  02*  and  NO*  molecules  at  high 
altitudes,  the  vertical  movement  of  the  luminosity  profiles  in 
response  to  the  passage  of  the  gravity  wave  is  much  smaller 
than  that  of  the  electron  density  profile.  The  amplitude  of  os- 
cillation of  the  airglow,  on  the  other  hand,  is  greater  than  that 
of  the  total  electron  density,  which  varies  monotonically  with 
time.  Owing  to  the  presence  of  a magnetic  field  as  well  as  the 
shape  of  the  electron  density  profile,  phase  differences  in 
behavior  occur  between  different  heights.  This  effect,  com 
bitted  with  the  height-differentiating  effect  of  quenching, 
produces  phase  differences  also  between  airglow  and 
ionospheric  parameters.  A combination  of  airglow  and 
ionospheric  nteasurmenls  can  be  used  to  differentiate  between 
gravity  waves  and  the  effects  of  incident  particle  flux. 

The  theory  has  been  tested  by  assuming  a gravity  wave  effect 
for  some  published  data  on  simultaneous  measurements  of  air- 
glow and  ionospheric  parameters.  The  fit  was  satisfactory. 

This  study  represents  the  first  theoretical  treatment,  as  far  as 
we  are  aware,  of  the  effects  of  gravity  waves  on  airglow  We 
have  tried  here  to  establish  the  basic  characteristics  to  be  ex 
peeled  for  the  luminosity  and  electron  density  profiles  undci 
the  influence  of  a gravity  w ave  and  to  show  how  these  could  be 
used  in  some  cases  to  provide  information  on  processes.  Ex- 
pansion of  the  techniques  should  be  useful  as  a tool  for  the 
study  of  both  gravity  waves  and  upper-atmosphere  parameters 
and  processes. 

APPt-.NDIX 

A good  analytical  approximation  for  the  decay  constant 
expressed  as  a function  of  the  scale  height  factor  p is  given  by 

= (B„Pn)wt 
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A derivation  for  this  is  given  by  Tuan  |l%8).  Another 
totally  different ‘derivation  has  been  given  by  Porter  and  hum 
|I974).  In  this  expression,  X„ decreases  monotonically  with  in- 
creasing />,  where  I < P <,  2.  Thus  the  decay  of  the  F layer  is  in 
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general  faster  when  the  gases  are  perfectly  mixed  (P  = I). 
The  expression  B(q)  in  (2)  is  given  by 


B(q)  = 


Aa“  sin  d cos  0I'[(3/2)  — g] 


As  was  already  mentioned,  q varies  between  0 and  'h.  If  q is 
set  to  zero,  the  velocity  amplitude  does  not  vary  with  height. 
Otherwise,  it  increases  exponentially.  We  see  that  B(q) 
vanishes  at  the  magnetic  pole  and  equator,  where  we  do  not 
expect  a horizontally  ducted  gravity  wave  to  produce  os- 
cillations in  airglow. 

The  functions  Pn(/>,  u.  t)  and  R„(p,  u>,  I)  are  defined  [ Porter 
and  Tuan.  1974]  by 

P»  = c(q,  t0,  t0)  J'  exp  — J 

•exp  — [X0(/'  — t0)  — B(q)  sin  u(t'  - /„)]  dt’  (A3) 

R„  = c(q,  t0,  i0)  J exp  — J 

•exp  -[X„(r'  - r„)  - B(q)  sin  oi(t’  - ?„)] 

• cos  a >(t'  — t„)  dt'  (A4) 
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k = H/2{B0D„?n 

Although  both  P„  and  R„  are  oscillating  functions  of  time  with 
frequency  w,  R„  may  change  sign,  but  P„  remains  always 
positive.  It  is  easy  to  see  that  as  w 0,  both  P„  and  /?„  become 
monotonic  functions  of  time,  and  all  oscillations  vanish  from 
the  expression  given  by  (2). 

The  functions  F„(p)  and  G„(q)  are  given  by 

. a"T(/»  - p)nn  + h)  T(n  — 1)T(1  + }) 

FM  = W rc=l> 


- . . I'(<7  + n - l)T[(3/2)  - q]  1(1  - q) 

Gn(q)  = — rr r — 


The  derivation  of  these  functions  is  given  by  Tuan  [1968] 
and  Porter  and  Tuan  [1974],  They  have  the  following  impor- 
tant properties: 

F„(P  = 1)  = 0 V n 

F.(p  = 2)  = 0 Vn  > 3 

GJq  = 0)  = 0 V/r  > 2 

Hence  the  infinite  sum  on  the  right-hand  side  of  (2)  will 

reduce  to  only  a few  terms  for  a gas  that  is  either  perfectly 
mixed  or  in  diffusive  equilibrium  and  for  a constant  velocity 
amplitude  (q  = 0). 
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